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ABSTRACT 


This  thesis  develops  a  means  by  which  a  closed  loop  frequency 
response  may  be  calculated  by  digital  computer,   given  only  the  circuit 
configuration  and  component  transfer  functions.     Using  these  data,    a 
program  is  developed  to  perform  the  necessary  logic  for  circuit  analy~ 
and  obtain  a  performance  function  and  frequency  response  over  any  de- 
sired frequency  range,   for  an  input  and  output  at  any  specified  circuit 
location. 

The  method  is  demonstrated  by  obtaining  the  performance  func- 
tion and  frequency  response  of  a  complicated  closed  loop  system. 
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OBJECT 


The  object  of  this  work  is  to  determine  the  feasibility  and 
devise  a  method  for  numerical  computation  of  closed  loop  fre- 
quency response  and  to  apply  the  method  to  an  illustrative  pro- 
blem. 
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CHAPTER  1 
INTRODUCTION 


In  the  design  and  analysis  of  feedback  control  systems  it  is 
often  necessary  or  desirable  to  know  the  steady  state  sinusoidal 
frequency  response  of  the  system.     Various  methods  for  obtaining 
this  frequency  response  are  available  to  the  engineer,    the  most 
often  used,   perhaps,   being  those  based  on  the  Bode  Diagram  and 
Nichols  Chart.     For  complicated  high  order  systems,    these  graph- 
ical techniques  are  apt  to  become  tedious  and  time-consuming;   even 
more  disadvantageous  is  that  the  time  must  be  spent  by  the  engineer, 
or  at  best,   an  assistant  with  knowledge  of  the  engineering  principles 
involved.     There  is  also  to  be  considered  the  human  factor  and  its 
accompanying  probability  of  error. 

It  would  be  advantageous  then  to  have  available  a  method  re- 
quiring a  minimum  application  by  the  engineer,    one  that  is  relatively 
fast,   free  from  the  human  factor  and,    after  the  minimum  initial  ap- 
plication,   could  be  handled  by  less-technically-qualified  personnel. 
With  the  advent  of  high  speed,   large  data  capacity  computers,    it  is 
reasonable  to  assume  that  such  a  method  could  be  made  available. 
It  is  toward  this  end  that  the  method  to  be  presented  in  this  paper  is 
directed. 

In  general,    the  minimum  input  information  required  from  the 
design  engineer  is  considered  to  be: 

1.  the  transfer  functions  of  the  individual  components  of 
the  system 

2.  a  means  of  describing  how  the  components  are  connected 

3.  the  frequency  range  and  desired  increment  of  frequency 
over  this  range  for  which  the  response  is  required. 

Item  1  is  assumed  known  to  the  engineer  and  need  only  be  entered  in 
the  form  described  in  Chapter  2.     Item  2  is  desired  to  be  as  simple  as 
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possible  but  still  be  of  such  a  nature  that  a  positive  check  may  be  made 
by  the  engineer  to  ensure  that  the  information  is  entered  correctly. 
This  is  described  in  detail  in  Chapter  2.     Item  3  consists  only  of  a 
statement  of  the  values  noted. 

It  is  noted  here  that  no  new  mathematical  concepts  are  set  forth, 
in  fact,    the  simplest  of  arithmetic  operations  are  involved.     These  basic 
fundamentals  are  the  most  readily  adapted  to  digital  computation. 
Graphical  solutions  such  as  the  Bode  Diagram  were  designed  not  be- 
cause the  fundamental  mathematics  were  unknown,  but  because  these 
simple  arithmetic  operations  were  so  laborious  as  to  be  prohibitive. 
The  modern  computer  performs  these  operations  at  an  extremely 
rapid  rate  and  permits  programming  the  heretofore  prohibitive  methods. 

Chapter  2  contains  the  physical  interpretation,   with  appropriate 
block  diagrams  and  flow  charts,    of  how  the  program  works;   i.e.  ,    how 
the  information  is  read  into  the  computer  and  the  general  steps  neces- 
sary to  compute  the  frequency  response.     The  physical  reasoning  and 
diagrams  are  intended  to  be  general  in  that  they  describe  the  necessary 
programming  steps  and  hence,   are  applicable  to  any  digital  computer. 
At  the  same  time,    they  are    in  great  enough  detail  to  permit  the  pro- 
gram to  be  adapted  to  any  computer  with  a  minimum  amount  of  additional 
labor. 

Chapter  3  and  appendices  thereto  contains  the  detailed  program 
as  written  for  the  IBM  650  Computer  and  its  associated  equipment  RAMAC, 
The  "MAC"  language  in  use  at  the  Instrumentation  Laboratory  at  MIT  is 
utilized  because  of  the  relative  ease  of  programming. 

Chapter  4  contains  a  brief  summary  of  the  salient  features  of  the 
program  and  its  application  to  some  sample  systems. 

The  primary  objective  of  this  thesis  is  to  set  forth  a  general  pro- 
gram,  stressing  simplicity  of  use  and  not  necessarily  to  derive  the  opti- 
mum program  with  respect  to  minimum  time  and  compactness.     For  this 
reason,   it  is  realized  that  the  program  of  Chapter  3  may  not  be  the  ulti- 
mate in  computer  programming  but  is  general  enough  to  be  adapted  to 
any  modern  digital  computer. 
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CHAPTER  2 
GENERAL  DESCRIPTION  OF  THE  SYSTEM 


2.1.  General 

The  program  that  is  used  to  determine  a  general  form  of  the 
system  performance  function  and  from  this  compute  the  frequency  res- 
ponse of  a  linear,    time  invariant  system  is  described  in  this  chapter. 
The  description  is  meant  to  be  quite  general  in  that  is  applies  to  pro- 
gramming any  digital  computer.     If  not  general,    the  program  would 
be  of  little  value  for  any  computer  except  the  specific  one  used  here. 
Even  though  general,   a  certain  amount  of  detail  is  necessary  to  en- 
able a  prospective  user  to  program  the  system  without  an  unduly  great 
amount  of  effort.     The  description  in  this  chapter  is  taken  by  steps, 
beginning  with  a  broad  over-all  aspect  and  block  diagram  flow  chart 
and  working  down  to  the  details.     This  gives  the  reader  an  opportunity 
to  decide  just  how  far  he  wants  to  go  and  permits  comprehension  at 
any  desired  level  of  detail. 

2.2.  Overall  Description 

As  stated  in  Chapter  1,    there  are  numerous  methods  in  the  litera- 
ture and  in  current  usage  to  calculate  the  frequency  response.     Some 
of  these,   particularly  the  graphical  techniques,    are  not  readily  adapt- 
able to  digital  computation.     We  may  note,   however,    that  in  obtaining 
the  system  performance  function  by  any  method,   we  get  at  some  stage 
of  the  calculation  an  equation  of  the  form: 


G„  G,    G^  +  GH  G     Gf       +  . 

a     b      c         d     e      t 

1  +  G^Gg  +  G2G3G4     +  . 
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where  the  G' s  represent  the  transfer  functions  of  the  system  components. 
A  more  careful  study  reveals  that  the  components  of  the  denominator 
terms  are  those  components  making  up  a  loop  or  combinations  of  loops 
in  the  physical  system.     Although  not  quite  so  obvious,    it  turns  out  that 
some  of  the  terms  of  the  numerator  may  be  recognized  as  loops  and  some 
as  direct  paths  from  the  input  to  the  output  of  the  system.     Therefore,   it 
would  seem  that  if  we  could  determine  the  general  scheme  of  how  these 
loops  and  paths  combine  in  the  equation,    then  it  might  be  possible  to 
readily  adapt  this  to  digital  computation.     Mason  (Ref.  5  )  obtains  the 
general  format  of  how  they  combine  and  arrives  at  a  general  equation 
representing  the  performance  function  for  an  input  at  i  and  an  output 
at  j  as: 

_J.    (S)  =  -a l i 


9, 


1  +  f  Lr  +  f  n2  Lr  +  f  H3  Lr  +  f  H  Lr  +  •    •    ■ 


where: 

P,  is  the  product  of  transfer  functions  of  all  components  in 

a  direct  path  from  i  to  j.     This  product  is  defined  as  the 

"Path  function"  of  path  "k". 

L  is  the  product  of  transfer  functions  of  all  components  of 

a  loop  and  is  defined  as  the  "loop  function"  of  loop  "r". 

nrt  L       is  products  taken  2,    3,    ...   at  a  time  of  loop  functions 
2      r  r 

of  all  loops  not  "touching"  one  another. 
L   ,  is  the  loop  function  of  a  loop  not  "touching"  path  "k" 

n    L        is  the  products  taken  2,    3,    ...   at  a  time  of  the  loop 
2     rk 

functions  of  all  loops  neither  "touching"  one  another  nor 

touching  path  "k" 
The  only  rule  involved  is  that  of  signs.     This  may  be  stated  as:   "In  trav- 
ersing a  loop,    if  the  number  of  (-)  signs  encountered  is  odd,    the  sign  of 
the  term  is  (+).     If  the  number  is  even  (including  zero),    then  the  sign  of 
the  term  is  (-)".     This  is  readily  understood  from  basic  feedback  system 
theory  where  we  are  familiar  with  the  result  that  a  negative  feedback 
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function  leads  to  a  characteristic  equation  1  +  the  loop  function.  In 
traversing  a  "direct  path"  if  the  number  of  (-)  signs  encountered  is 
even  (including  zero)  the  sign  of  the  term  is  (+);  if  odd,  the  sign  of 
the  term  is  (-). 

To  demonstrate  the  use  of  the  basic  equation  and  the  "touching 
vs.   non-touching"  character  of  loops  and  paths  consider  the  following 
system  where  G     is  the  transfer  function  of  block  (N) 


From  the  block  diagram,    it  is  seen  that  there  are: 
1)    4  loops  with  loop  functions 


Ll  =  GlG2G3G4G6 
L2  =  G5G4G6 
L3  =  G2G7 
L4  =  "  G3G4G8 


S  L  =  G,G„G0G„G„  +  GrG.G,,  +  G^G,,  -  G0G„G 


'1^2  3  4~6 


5~4~6 


'2~7 


'3~4~8 


2)    2  direct  paths  with  path 
functions 


Pl  =  GlG2G3G4 
P2  =  G5G4 
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The  "touching"  of  2  loops  can  also  be  seen  by  inspection  but  if  in  doubt 
the  question  might  be  asked  "can  a  block  be  inserted  in  one  of  the  loops 
without  changing  the  loop  function  of  the  other?"    If  so,    then  they  do  not 
touch.     In  the  block  diagram  above: 

loop  (1)  touches  the  remaining  three 

loop  (2)  does  not  touch  loop  (3)  but  does  touch  (1)  and  (4) 

loop  (3)  and  (4)  do  not  touch  one  another 


then; 


S   n2  Lr  =  L2  L3  +  L3  L4  -  G5G4G6G2G?  -  G^G^Gg 

Path  (1)  is  seen  to  touch  all  loops:  ZL    .   =0 

Path  (2)  touches  all  except  loop  (3):        2  L   2  =  G2G? 

The  performance  function  then  is: 

qout(S)_ G1G2G3G4  +  G5G4  +  G5G4G2G7 

1  +G1G2G3G4G6+G5G4G6+G2G7  -  G3G4Gg+G5G4G6G2G7  -  G^^G^g 

To  obtain  the  frequency  response,    two  courses  are  now  available. 
The  component  transfer  functions  might  now  be  entered  and  a  program 
written  to  expand  the  performance  function  above  into  the  form 

ASn  +  BSn_1  +  CSn"2  +.,. 


_  cm  cm-l  cm-2 , 

aS      +bS  +cS         jr  .  .  . 


Then  let  S  =  j  U  and  compute  the  desired  response.  This  method  would 
require  extensive  working  and  storage  space  and  might  well  exceed  the 
capacity  of  some  computers.     The  other  course  available  is  to  enter 
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the  transfer  functions  of  the  components,   let  S  =  j   u    and  compute  a 
magnitude  ratio  and  phase  angle  for  each  component  for  a  particular 

value  of  to.     Then,    substitute  these  values  into  the    — r-     equation, 

q  in       n 

perform  the  necessary  mathematical  operations  and  obtain  the  res- 
ponse for  the  particular  value  of  u.     The  program  described  in 
Chapter  3,   uses  the  latter  method. 

The  form  in  which  the  final  answer  is  received  will  vary  with 
the  computer;   perhaps  punched  cards,    or  if  a  plotter  is  available, 
then  the  final  result  may  be  read  directly  from  a  plot.     A  flow  chart 
of  this  broad  description  is  contained  in  Fig.    2.  1. 


2.3.     Input  Data 

It  is  desired  that  this  be  in  as  simple  a  form  as  possible  and 
provide  some  sort  of  a  check  to  the  programmer  that  the  informa- 
tion (data)  has  been  correctly  entered.     To  fulfill  both  of  these  re- 
quirements,   it  is  considered  best  to  use  an  arbitrary  numbering 
system  of  all  components  and  juction  points  (summers).     Then  the 
inputs  to  a  numbered  unit  are  the  numbers  of  the  unit  (s)  from  which 
they  came.     All  that  is  necessary  to  be  entered  then  is  the  number  of 
the  unit  and  the  number  (s)  of  the  input.     To  permit  operation  by  the 
computer,    these  data  are  entered  in  an  array.     For  example,    con- 
sider the  simple  system  illustrated  here  in  block  diagram  form: 

4 

out 


If  computer  space  is  critical  or  if  the  system  is  extremely  large,    the 
blocks  2  and  3  may  be  just  as  well  combined  into  a  single  unit.     If  con- 
sidered as  two  separate  units,    the  input  array  would  be: 

Unit  Inputs 

2  4 

3  2 

4  1-3    (the  1  is  unnecessary  but  used  here  to 

illustrate  the  form  for  multiple  inputs). 
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JL 


Read  Input  Data 


Compute  Performance 
Function  of  System 


Set  First 
Frequency 


Compute  Magnitude  Ratio 
And  Phase  Angle  For  One 
Frequency 


f  No     ,<- 


Is  This  The 
Last  Frequency? 


->f    Yes  j 


2k. 


Increase  Frequency 
by  Increment 


Fig.    2.  1.     Flow  Chart  of  Broad 

Aspect  of  the  Program 
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Note  the  (-)  indicating  this  is  to  be  subtracted  in  the  summer;    i.e.   neg- 
ative feedback. 

The  only  other  input  data  required  are:     input  and  output  location  num- 
bers,   in  this  case  4,  3;   number  of  components  and  summers  and  max- 
imum number  of  inputs  to  any  unit.     The  reason  for  the  first  require- 
ment is  obvious;    the  last  three  are  used  to  reduce  computation  time 
and  conserve  computer  storage  space  by  establishing  the  exact  size  of 
the  array. 

For  reasons  to  be  explained  later,   it  is  necessary  to  number  the 
summers  also  and  is  convenient  to  assign  larger  numbers  to  these 
than  the  highest  numbered  component  and  list  them  last  in  the  array. 
The  particular  order  in  which  either  components  or  summers  are 
numbered  is  otherwise  arbitrary. 

2.4.     Computation  of  Performance  Function 

Using  the  method  described  in  paragraph  2.2.  ,   it  is  seen  that 
what  is  actually  required  is  a  method  of  tracing  through  the  numerical 
input  data  to  determine  the  existent  loops  and  direct  paths  of  the  sys- 
tem.    These  are  found  and  recorded  as  groups  of  numbers;   hence  to 
determine  the  "non-touching"  character,    it  is  necessary  only  to  check 
numbers  of  each  group  against  the  other  groups.     Since  it  is  possible 
for  paths  or  loops  to  "touch"  without  having  common  components,    the 
necessity  for  numbering  the  summers  is  now  seen.     It  is  not  possible 
for  loops  to  touch  without  having  at  least  one  common  junction  point 
(summer)  and  vice  versa.     The  broad  scheme  of  this  operation  is 
described  in  the  flow  chart  of  Fig.   2.2. 

After  determining  the  existent  loops,   paths,    non-touching  loops, 
etc.  ,    the  groups  of  numbers  are  compiled  to  form  the  general  expres- 
sion for  the  performance  function.     Since  the  requirement  might  exist 
for  knowing  the  performance  function,   a  subroutine  is  included  in  the 
main  program  for  compiling  and  punching  out  this  function,   before 
commencing  calculation  of  frequency  response. 

The  methods  usted  to  determine  the  paths  and  loops  are  perhaps 
the  heart  of  the  entire  program.     The  two  methods  are  similar  in 
general  hence  only  the  loop  routine  will  be  explained  in  detail  and  the 
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Read  Input  Data 

v 

Determine  all  Loops  (L    ) 
of  the  System 

V 

Determine  all 
Direct  Paths  (Pk) 

Determine  all  Products 
of  non-touching  loops  (IlL   ) 


Determine  all  loops  not 
touching  each  Path  (L   ,) 


^U- 


Determine  all  Products  of 
non-touching  loops  that  do 
not  touch  each  Path  (nL   ,  ) 


Combine  as  Indicated  in  the 
Basic  Equation  if  Perfor- 
mance Function  desired. 
Otherwise  store  in  Form  de- 
termined above. 


Fig.   2.2.     Calculation  Of  The  Performance  Function 
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path  routine  in  general  except  where  it  differs.     The  first  step  of  the 
loop  routine  is  to  choose  a  starting  point  in  the  array.     Since  any  loop 
of  the  system  must  contain  a  junction  point  (summer),    then  we  choose 
as  the  starting  point  the  first  summer  listed  in  the  input  data  array 
(not  necessarily  the  first  one  in  the  physical  system).     We  could  start 
with  a  component  but  can  save  time  by  considering  only  the  summers. 
Hence,   it  is  convenient  to  list  the  components  first  then  the  summers. 
This  order  might  well  be  reversed  but  the  program  of  Chapter  3  is 
written  for  the  order  noted.     The  first  summer  is  located  by  knowing 
the  number  of  components  other  that  summers.     To  enable  computa- 
tion of  the  frequency  response  or  in  punching  the  performance  function, 
the  summer  numbers  are  not  desired.     This  is  obviated  by  assigning 
higher  numbers  to  the  summers  than  to  the  components,    and  discard- 
ing all  numbers  greater  than  that  of  the  highest  numbered  component. 
The  computer,   after  noting  the  starting  number,   then  scans  the  input 
array  by  rows  until  this  number  is  found  as  an  input  to  another  num- 
bered unit.     Then  the  new  number  is  used  as  the  "search  number" 
until  it  is  found  as  the  input  to  another  numbered  unit,   and  so  on  until 
we  have  a  complete  loop.     It  is  seen  now  that  difficulty  might  be  en- 
countered due  to  the  output  of  a  unit  going  to  more  than  one  other  unit. 
This  must  be  considered  in  the  program.     The  way  this  is  done  is  to 
record  the  first  output  then  continue  the  search  to  see  if  there  are  any 
more.     If  more  are  found  they  too  are  recorded  at  this  time.     It  is 
impractical  to  find  one  complete  loop  then  go  back  to  the  starting  point 
and  search  for  more  loops.     Hence,  all  points  in  all  branches  are  re- 
corded at  this  time.     This  is  more  clearly  described  by  the  flow  chart 
of  Fig.   2.3.     After  recording  (in  the  computer)  all  the  loops  contain- 
ing the  first  summer,   the  program  continues  to  the  second  and  re- 
peats the  search  and  record  process,   ensuring  that  loops  already  re- 
corded are  not  repeated,   and  so  on  until  all  loops  are  found.     It  is 
realized,   of  course,   that  checks  must  be  made  for  such  things  as 
traversing  a  loop  only  once  and  recording  only  complete  loops.     This 
is  also  brought  out  more  clearly  in  Fig.   2.3. 
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Fig.   2.3.     Detailed  Flow  Chart  of  Loop  Routine 
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The  path  routine  differs  from  the  loop  routine  primarily  in  that 
the  starting  point  is  the  number  inserted  in  the  input  data  as  the  sys- 
tem input  for  which  the  response  is  desired,   and  also  the  end  of  the 
path  is  the  listed  output  number.     The  similarities  and  differences 
are  seen  more  clearly  by  comparing  the  flow  chart  for  the  path  rou- 
tine (Fig.    2.  4.)  with  Fig.   2.3. 

The  importance  of  these  two  subroutines  and  the  program 
sophistication  necessary  to  accomplish  the  desired  result  cannot 
be  emphasized  too  strongly.     The  flow  charts  of  Figs.   2.3.   and  2.4. 
clarify  the  procedure  satisfactorily  for  understanding  the  general 
method  employed.     However,    they  are  not  considered  to  be  in  great 
enough  detail  to  permit  writing  a  similar  program  for  another  com- 
puter.    For  this  reason,    the  method  of  this  subroutine  is  explained 
in  great  detail  in  Chapter  3. 

Having  recorded  all  the  existent  loops  and  direct  paths  of  the 
system,    it  remains  to  determine  the  remaining  factors  of  the  general 
equation.     The  idea  is  very  simple;   we  merely  check  numbers  of 
groups  already  recorded  against  appropriate  other  groups.     If  no 
identical  numbers  are  found  then  the  loops,   etc.  ,   do  not  touch.    Since 
this  procedure  is  fairly  readily  understood,    a  general  flow  chart 
(Fig.    2.5.)  is  considered  adequate. 

To  obtain  the  performance  function,   we  now  combine  the  ap- 
propriate previously  stored  numbers  and  punch  out  the  function  term 
by  term.     Fig.    2.  6.  ,    is  a  general  flow  chart  of  this  routine.     It 
should  be  pointed  out  that  since  it  is  foreseeable  that  the  performance 
function  might  not  actually  be  desired,    then  writing  this  subroutine 
is  optional.     For  this  reason  a  "by-pass"  is  put  into  the  program  and 
effected  by  entering  a  value  of  zero  for  a  variable  read  in  with  the 
input  data.     If  the  performance  function  is  desired,    then  any  number 
other  than  zero  is  entered. 

2.5.     Component  Transfer  Functions  and  Frequency  Response 

Any  component  whose  input-output  performance  can  be  described 
by  a  linear  ordinary  differential  equation  with  real  constant  coefficients 
has  a  transfer  function  which  is  rational  in  the  Laplace  variable  S 
with  the  numerator  and  denominator  polynomials  expressible  as 
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Fig.    2.4.     Detailed  Flow  Chart  of  Path  Routine 
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2 
products  of  terms  of  the  form  AS     +  B  S  +  C  in  which  any  of  the 

coefficients  may  be  zero.     These  transfer  functions  must  now  be 

entered.     The  form  in  which  they  are  entered  is  discussed  in  detail 

with  the  program  of  Chapter  3  and  illustrated  in  Chapter  4. 

Next  a  magnitude  ratio. and  phase  angle  is  computed  for  this 
component  by  letting  S  =  j   w    and  substituting  the  particular  value 
of  w  for  which  we  are  computing  the  response.     The  required  bits 
of  information  regarding  frequency  are:     frequency  range  (upper 
and  lower  limit)  and  increment.     (See  Fig.    2.  7.  ). 

To  facilitate  ease  of  recording  the  required  information,   a 
form  is  included  with  the  program  details  of  Chapter  3  and  the 
sample  system  of  Chapter  4. 

2.  6.     System  Magnitude  Ratio  and  Phase  Angle 

Having  computed,   for  the  particular  value  of  w,    the  magni- 
tude and  phase  angle  of  each  component  the  next  and  final  step  is  to 
substitute  these  values  into  the  general  performance  function.     If 
the  entire  performance  function  has  been  retained  in  storage,    if 
and  when  punched,   it  could  well  be  used  to  compute  the  system  res- 
ponse;  but,    since  it  is  desired  to  have  this  performance  function 
calculation  optional,    then  the  routine  for  finding  system  response  is 
written  using  the  loops,   paths  etc.  ,   previously  computed.     This  is  a 
straight -forward  routine  easily  written  and  understood  once  it  is 
known  where  the  information  is  stored  from  previous  subroutines. 
Since  the  general  equation  from  paragraph  2.  2.  ,    contains  only  loops 
and  paths  and  combinations  of  these,    it  is  expedient  to  compute  and 
store  first  the  magnitude  and  phase  for  each  loop  and  path,    then  use 
these  stored  values  in  the  combinations.     The  performance  function 
is  in  general  a  summation  of  loops,    paths,   etc.  ,   hence  a  conversion 
from  polar  to  rectangular  form  is  also  required  and,   after  compil- 
ing,  a  conversion  back  to  polar  form.     For  any  greater  detail,    the 
reader  is  referred  to  the  program  detail  of  Chapter  3  and  the  sample 
problem  of  Chapter  4. 
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Fig.   2.7.     Flow  Diagram  of  Component 

Magnitude  and  Phase  Anjgle  Routine 
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2.7.     Frequency  Response  "Re-Run" 

After  obtaining  a  frequency  response  over  an  indicated  range 
with  a  specified  increment,    it  might  be  desired  to  change  the  range, 
increment  or  even  the  transfer  function  of  some  of  the  existent  com- 
ponents.    In  this  case,   it  would  not  be  desired  to  recompute  the  gen- 
eral performance  equation  of  paragraph  2.  2.  ,   hence,   provisions 
are  made  to  re-enter  the  program  with  these  different  values.     The 
method  will  vary  with  the  computer  used  anti  for  this  reason,   is  not 
detailed  here.     A  detailed  step  by  step  procedure  for  the  IBM  650 
is  contained  in  paragraph  3.10. 


-21- 


CHAPTER  3 
DETAILED  DESCRIPTION  OF  PROGRAM  WRITTEN  FOR  IBM- 650 


3.1.     General 

It  was  stated  in  Chapter  1,    that  the  MAC  language  (Ref.    6)  is 
used  in  this  program  for  simplicity  of  construction.     Very  briefly, 
MAC  may  be  described  as  follows:     A  plug  board  wiring  scheme  per- 
mits recognition  of  certain  special  plain  language  commands  such  as 
"PUNCH",    "GO  TO",    "IF",    "DO"  and  "DO  TO".     Arithmetic  symbols 
are  also  recognized;   examples  of  these  are  the  signs  for  "equals", 
"plus",    "minus",    "decimal  point"  and  "parenthesis".     Any  number  of 
parentheses  may  be  used  consistent  with  proper  arithmetical  usage. 

Variables  are  created  by  several  means,    the  most  common  one 
being  an  equation  with  the  variable  to  the  left  of  the  equal  sign,   e.g.  , 
A  =  2.  5  4-    B/2  is  interpreted  as;   "a  new  variable  (or  new  value  of  a 
variable  in  use)  is  equal  to  2.  5  plus  (B  (a  previously  defined  quantity) 
divided  by  2).     Upon  encountering  this  command  in  a  program,    the 
value  of  A  is  computed  and  stored  in  the  memory  and  may  be  used  at 
any  time  in  the  right  hand  side  of  an  equation. 

Routines  for  certain  other  functions,    such  as  trigonometric 
functions,    square  root,    absolute  value,    integer  are  built  into  the  basic 
storage  and  are  referred  to  upon  recognition  of  the  functions:   SIN, 
COS,    TAN,   ARCSIN,   SQRT,   ABS,   INTEGER  and  others  for  which  the 
reader  is  referred  to  Ref.    6  of  the  bibliography.     Any  number  of  sub- 
scripts,  either  numbers  or  letters  for  which  a  value  has  been  compu- 
ted,   may  also  be  used.     It  is  noted  that  subscripted  variables  are  treat- 
ed in  the  same  manner  as  non-subscripted  ones  and  new  variables  may 
be  introduced  accordingly,     If  the  subscript  is  itself  a  variable,    its 
value  must  be  defined.     The  subscript  is  then  an  "index"  and  must  be 
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so  indicated  on  the  first  card  of  the  program.     Variables  are  limited 
to  letters  of  the  alphabet,   numbers,   and  combinations  of  either  or 
both.     The  physical  limits  imposed  by  MAC  is  a  total  of  800  variables. 

The  individual  lines  and  cards  of  the  MAC  program  are  number- 
ed in  sequence,    using  the  prefixes  E,   M,   S  to  indicate  exponent,    main, 
and  subscript  cards  respectively,     Equation  numbers  when  required 
appear  to  the  right  of  the  MAC  statements  and  equations,    and  are  used 
as  a  means  of  transferring  control  to  selected  points  in  the  program 
or  to  establish  loops. 

After  all  MAC  program  cards  are  read  into  the  computer,   an 
internal  assembly  routine  converts  this  program  into  the  basic  language 
of  the  computer,   begins  at  a  designated  point,   and  proceeds  with  the 
computation.     This  "new  program"  may  be  punched  on  cards  as  it  is 
being  written,    the  cards  then  being  referred  to  as  a  "load  deck";   or 
if  no  future  need  is  anticipated,   the  load  deck  is  not  punched.     Use  of 
this  load  deck  in  lieu  of  MAC  cards  for  future  operations  reduces  com- 
putation time  considerably. 

3.2.     Major  Subroutines  of  the  Program 

Reference  to  the  basic  equation  for  determining  the  performance 
function,   Paragraph  2.2.,   shows  that  the  component  numbers  are  re- 
quired for: 

a.  All  closed  loops  of  the  system  (L   ) 

b.  All  direct  paths  from  input  to  output  (P/) 

c.  All  combinations  of  2,    3,    4  etc.   loops  not  touching  one 

another  (L.     products) 
r 

d.  All  loops  not  touching  each  path  (L.   ,  ) 

e.  All  combinations  of  2,    3,    4  etc.   loops  not  touching  each 
other  and  not  touching  each  path  (L   ,    products) 

In  addition  to  the  above,   other  major  subroutines  are: 

f.  Calculation  and  punching  of  the  performance  function 

g.  Calculation  of  frequency  response  (for  each  w)  of: 

1.  All  components 

2.  All  loop  functions,    all  path  functions,    certain  pro- 
ducts of  loop  and  path  functions 

3.  Denominator  of  the  system  response 

4.  Numerator  of  the  system  response 
5„     The  system 
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Flow  charts  for  each  of  these  have  been  included  in  the  general  dis- 
cussion of  Chapter  2.     The  complete  program  is  appended  as  Appen- 
dix B  and  is  discussed  in  detail  in  the  following  paragraphs.     Each     . 
major  subroutine  will  be  explained  separately,   with  reference  to  the 
appropriate  flow  chart  and  sections  of  the  program.     The  numbers  in 
parentheses  following  the  subroutine  paragraph  title  indicate  the  appro- 
priate location  in  the  program. 

3.3.     Input  Data 

Before  discussing  the  program,  it  is  essential  that  the  "Input 
Data  Array"  be  understood.  This  is  most  readily  accomplished  by 
means  of  a  block  diagram  of  a  hypothetical  system: 


Fig.    3.1. 

The  numbering  of  components  and  summers  is  shown  to  be  purely 
arbitrary  except  that  summers  are  assigned  higher  numbers.     Data 
Cards  for  MAC  contain  six  or  less  numbers  and  are  written  in  floating 
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point  form  (Ref.   7).     The  first  card  contains  in  the  order  listed: 
Maximum  no.   of  inputs  to  any  unit 

(excluding  the  system  input)  2 

Number  of  component  boxes  5 

Number  of  summers  3 

Input  unit  number  7 

Output  unit  number  2 

Is  performance  function  desired  0  for  No. 

1  for  Yes 

On  the  succeeding  cards,  the  inputs  to  each  unit  are  listed,  using  a 
separate  card  for  each  unit.  Again  the  order  of  listing  is  arbitrary 
except  that  the  summers  should  follow  the  component  boxes: 

s  unit 

(0-3)  (numbers  in  parentheses 

(4-7)  are  not  entered,   but  are 

(8-11)  used  for  explanation  fol- 
(12-15)  lowing) 
(16-19) 
(20-23) 
(24-27) 
(28-31) 
(32-35) 

The  zeros  are  required  since  with  MAC  a  "Read"  command  is 
used  to  enter  data  and  this  Read  command  must  specify  what  values  and 
how  many  are  to  be  read.     The  program  for  which  this  Chapter  is  writ- 
ten uses  the  "maximum  number  of  inputs  to  any  unit"  as  the  basis  for 
specifying  the  width  of  the  array  and  hence,    the  number  of  values  to  be 
read.     The  zeros  are  used  as  "stops"  in  the  search  routine.   Note,    the 
(-)  signs  indicating  negative  feedback  into  unit  7  and  8. 

For  any  system,   a  rectangular  array  is  obtained,    the  width  of 
which  is  the  "maximum  number  of  inputs  to  any  unit"  +  2.     The  length 
of  the  array  is  "the  number  of  component  boxes  +  the  number  of  sum- 
mers +  1  (the  last  row  is  all  zeros). 


Unit  No. 

1 

nputs 

to 

1 

8 

6 
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6 

0 
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7 

0 
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Each  position  in  the  array  is  identified  as  a  subscripted  variable 
and  the  values  read  from  the  data  cards  are  thereafter  recognized  as 
the  value  of  the  variable.     The  numbers  in  parenthesis  are  the  sub- 
scripts identifying  the  locations.     G  is  chosen  as  the  variable  to  be 
subscripted.     For  example,    the  value  of  G22  is  3. 

3.4.    L     Subroutine  (M  0075  -  M  0585) 
*r 

It  is  desired  to  determine  the  units  comprising  each  of  the  closed 
loops  of  the  system.     As  noted  in  Chapter  2,    this  is  a  very  important 
but  complicated  subroutine.     To  ensure  complete  understanding  requires 
that  the  program  be  explained  in  great  detail.     Since  this  is,   of  neces- 
sity,  a  lengthy  discussion  the  program  details  are  included  as  Appendix 
A  and  only  the  method  will  be  discussed  here.     Refer  to  the  flow  chart 
of  Fig.   2.  3.     The  basic  method  is  to  choose  as  a  starting  point  in  the 
array,    the  first  summer  listed.     Using  this  number  as  a  "search  num- 
ber" the  input  array  is  scanned,    starting  with  G,   until  the  number  is 
found  as  an  input  to  another  unit.     The  search  number  and  the  number 
of  the  output  unit  are  entered  temporarily  in  an  array  designated  by  the 
subscripted  variable  XX.     The  search  continues  for  more  outputs  of 
the  "search  number".     If  any  more  are  found,   a  new  line  is  formed  in 
the  XX-array  for  each.     After  all  outputs  are  recorded,    the  last  output 
is  chosen  as  a  new  search  number  and  the  same  procedure  carried  out 
again,   writing  a  new  line  in  the  XX-array  each  time  a  new  unit  is  found 
to  be  a  point  in  the  loop.     Although  not  too  obvious  at  this  time,    this  is 
necessary  to  ensure  the  same  path  is  not  being  traversed  around  a  loop, 
and  to  have  all  numbers  available  when  the  loop  is  found  to  have  been 
completed. 

To  illustrate  briefly,   the  first  summer  listed  in  the  array  of 
paragraph  3.3.,   is  (6)  which  goes  to  (2),   so  the  XX-array  is: 

line  (a)  6 

line  (b)  6    2 

No  more  outputs  of  6  are  found  so  using  (2)  as  a  search  number,    its 
output  is  found  to  go  to  (5)  and  negatively  to  (7).     These  are  recorded 
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in  the  XX  array  as  lines  (c)  and  (d): 


line  (c)  6      2      5 

line  (d)         -6      2      7 

(7)  is  the  next  search  number,   for  which  an  output  is  found  to  3  and  4. 
These  are  recorded  as  lines  (e)  and  (f).     Then  (5)  is  selected  as  a 
search  number,   the  output  of  which  feeds  negatively  to  (8),    (line  g) 

line  (e)  -6273 
line  (f)  -6274 
line  (g)         -6258 

The  routine  continues,    checking  each  new  number  to  see  if  it  forms 
a  completed  loop.     If  so,    the  appropriate  line  of  the  XX-array  is  re- 
corded permanently  as  a  loop,   for  which  the  subscripted  variable  ZZ 
is  used.     Continuing  the  search  process  yields  one  new  line  from  (g) 
and  one  from  (f)  but  when  (3)  from  line  (e)  is  checked,   it  is  found  to 
go  to  (6)  which  is  the  starting  point,   hence  a  complete  loop  has  been 
found  and  a  new  line  is  not  entered  in  the  XX-array,   but  is  recorded 
as  ZZ,   to  ZZ.  with  values  of  -6,   2,    7,    3  respectively.     ZZ5  is  then 
recorded  as  zero  as  a  separator  from  the  next  loop. 

After  each  number  is  located,   it  must  be  checked  to  see  that  it 
has  not  already  been  entered  as  a  point  in  the  line  of  the  XX-array. 
If  so,    this  represents  a  minor  loop  within  the  main  loop  that  is  being 
traced.     This  minor  loop  will  be  recorded  as  a  line  in  a  succeeding 
array.     Hence,   no  new  line  is  entered.     This  is  illustrated  as  follows 


line  (h) 

-6 

2 

5 

8 

1 

line  (i) 

-6 

2 

7 

4 

8 

line  (j)         -  6     2      7      4      8      1 

Now  use  (1)  from  line  (h)  as  a  search  number  and  find  it  goes  to  (6), 
a  complete  loop,   and  back  negatively  to  (7).     A  new  loop  ZZfi  to  ZZ..  „ 
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is  recorded  as  -6,    2,    5,    8,    1  and  a  new  line  (k)  is  recorded  as: 

line  (k)        +625817      (note  the  sign  change) 


line  (1) 

6 

2 

5 

8 

1 

7 

3 

line  (m) 

6 

2 

5 

8 

1 

7 

4 

Lines  (1)  and  (m)  are  recorded  from  line  (k)  and  we  move  back  to  the  (1) 
of  line  (j).     (1)  is  found  to  go  to  (6),    a  complete  loop,   hence,   is  record- 
ed in  the  ZZ  list  but  not  in  the  XX-array,    and  to  (7);  but  since  (7)  has 
already  been  recorded  in  this  line,   it  is  disregarded  and  nothing  is  re- 
corded in  the  XX-array.     The  process  continues  until  all  loops  of  which 
summer  (6)  is  a  part,   have  been  recorded  in  the  ZZ  list.     The  next 
summer  of  the  "input  data  array"  is  then  chosen  as  a  starting  point  and 
the  same  procedure  is  carried  out  as  for  the  first  summer.     It  is  realiz- 
ed now  that  loops  that  have  already  been  recorded  from  summer  (6)  must 
not  be  recorded  again,    so  appropriate  checks  must  be  made  to  prevent 
this.     When  all  summers  have  been  checked  for  new  loops,   all  of  the 
existent  loops  of  the  system  have  been  recorded  and  the  program  moves 
on  to  the  next  major  subroutine.     For  greater  detail  as  to  how  the  above 
methods  are  programmed,    refer  to  Appendix  A. 

3.  5.     Subroutine  for  Determining  Direct  Paths  (PK  routine;   MO590-1040) 

This  routine  differs  only  slightly  from  the  LR  routine  above,    and 
will  be  discussed  in  detail  only  where  major  differences  occur.     The 
search  of  the  G-array  is  the  same  as  above  except  the  idea  now  is  to  find 
all  components  in  a  path  between  two  specified  units.     Therefore,    the 
first  number  of  the  XX-array  is  the  designated  input  unit  and  the  array 
is  written  only  once. 

Recording  new  lines  in  the  XX-array  (MO790-0820)  and  the  checks 
made  prior  to  this  (MO735-0785)  are  the  same  as  before  except  here,    the 
check  for  a  complete  path  is  made  against  the  specified  output  and  the 
duplicate  point  check  includes  the  starting  number.     When  complete  paths 
are  found,    all  units  of  the  path  are  recorded  as  the  subscripted  variable 
YY  and  this  procedure  differs  from  the  LR  routine  only  in  that  the  last 
number  found  is  also  recorded  (MO970). 
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It  is  necessary  to  know  the  number  of  paths  just  as  it  was  neces- 
sary to  know  the  number  of  loops,   hence  a  tally  is  made  each  time  a 
new  path  is  recorded  (M1020  and  M1040).     Note  that  YQ  is  the  number 
of  paths  and  the  subscripted  variable  Y  also  is  used  in  the  XX-array. 
The  reasoning  here  is  that  the  variables  used  in  writing  the  XX-array 
are  not  needed  after  this  subroutine  is  completed  and  with  limited  com- 
puter space  it  is  best  to  use  variables  that  may  also  be  used  later. 
Throughout  the  program,   the  single  letter  variables  are  used  many 
times  but  space  is  reserved  for  a  designated  number  of  variables  with 
indexed  subscripts  (i.e.  ,    subscripts  that  are  also  variables)  and,   once 
used,    this  space  is  wasted  unless  the  subscripted  variables  are  chosen 
such  that  they  may  be  used  again.     Since,   as  will  be  seen  later,    there 
are  a  number  of  constants  associated  with  each  path,    (e.g.   the  number 
of  loops  not  touching  each  path),   it  is  convenient  to  use  an  indexed- 
subscripted  variable  to  denote  these  constants.     The  subscripted  vari- 
able Y  is  used  for  this  and  hence  that  space  used  by  the  Y?  s  in  writing 
the  XX-array  is  not  wasted.     The  only  possibility  of  an  overlap  is  Yfi. 
Therefore,    the  tally  of  paths  is  kept  as  a  single  letter  variable  (M1020) 
and  upon  completion  Y~  is  set  equal  to  this  tally. 

3.6.     Products  of  Loop  Functions  Subroutine  (M1050-1795) 

From  the  basic  performance  function  equation  of  Chapter  2,   it  is 
seen  that  both  the  numerator  and  denominator  are  made  up  of  a  series 
of  terms,    the  number  of  which  is  determined  by  the  complexity  of  the 
physical  system.     It  is  impractical  to  program  an  unlimited  number  of 
product  terms  and  the  question  quite  naturally  arises  as  to  how  many 
should  the  program  be  capable  of  handling.     Intuitively,    it  is  seen  that 
the  degree  of  complexity  of  the  physical  system  also  determines  the 
amount  of  computer  space  required  for  the  "input  data  array",   the 
XX-array,    the  ZZ  and  YY  lists,   etc.     There  is  therefore,   a  physical 
limit  imposed  on  any  program  written  in  MAC  without  resorting  to 
other  means  for  obtaining  computer  space.     Consideration  of  these 
limiting  factors  leads  to  an  estimation  of  the  maximum  number  of  pro- 
duct; terms  that  should  be  programmed.     The  program  appended  to  this 
report  is  written  for  a  maximum  of  "products  of  loop  functions  taken 
four  at  a  time  of  all  loops  not  touching  one  another". 
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The  method  used  to  determine  the  touching  vs.    non-touching 
character  of  the  loops  is  to  check  each  number  (component  or  summer) 
recorded  in  a  loop  against  each  number  of  every  other  loop.     If  no 
identical  numbers  are  found  in  the  two  loops,    then  they  do  not  "touch" 
and  are  to  be  so  recorded.     The  details  of  the  program  are  most 
readily  explained  by  means  of  an  example. 

Consider  a  system  containing  four  loops.  These  have  been  re- 
corded in  the  ZZ  list  with  a  zero  separating  the  numbers  of  one 
loop  from  the  next: 


(1) 

(2) 

(3) 

(4) 

zzt 

ZZ 

a  + 

2 

ZZb  +  2 

ZZc  +  2 

zz2 

■ 
» 

• 

• 
• 

• 

• 

• 

• 

ZZ 

a 

ZZ, 

b 

ZZ 

c 

zzd 

ZZ 

a  + 

1 

0 

ZZb  + 

1  = 

0 

ZZ    1  = 

c  +  1 

0 

zzd  + 1 = 

0 

Check  loop  (1)  against  loop  (2)  by  first  checking  ZZ1   against  ZZ         9 

to  ZZ,     ,  1  . 
b  +  1 

For  convenience  of  locating  loops  by  number  at  a  later  time, 

the  first  ZZ  number  of  each  loop  is  recorded  as  the  check  proceeds. 

The  subscripted  variable  G  is  used  for  this  since  at  the  present  time, 

the  space  reserved  for  the  G' s  is  unused.     Since  a  +  2  is  not  yet  known, 

it  is  found  by  checking  the  ZZ  list  until  a  zero  is  found  (Ml  135-1 140). 

The  next  ZZ  number  after  the  zero  is  ZZ      ,   n  and  G0  is  recorded  as 

a  +  2  2 

this  value  (Ml  160).     Now  the  check  consists  of:     ZZ..   (first  check 
point)  against  all  numbers  of  loop  2.     When  ZZ,     .   1    (=0)  is  found, 
pick  up  ZZ0  as  the  "check  point"  and  start  again  with  ZZ      ,   0  as  the 
"check  number".     The  procedure  is  continued  until  ZZ         1    (=0)  is 
found  as  the  check  point  or  until  the  check  number  and  check  point  are 
identical.     Since  algebraic  signs  are  of  no  significance,   absolute  values 
are  used  in  the  check. 

If  the  check  point  becomes  zero,   loops  (1)  and  (2)  are  recorded 
as  a  pair  of  non-touching  loops,   and  loop  (1)  is  then  checked  against 
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loop  (3)  and  so  on.     When  the  loop  number  of  the  loop  to  be  checked 
exceeds  Z0  (number  of  loops)  (Ml  150)  the  procedure  is  to  check  loop 
(2)  against  the  remaining  loops  in  the  same  manner  as  before,   except 
now  the  first  ZZ  number  of  each  loop  is  known  and  used  as  the  first 
"check  number".     The  procedure  continues  until  loop  number  (ZQ  -  1) 
is  checked  against  loop  number  (ZQ). 

The  non-touching  loop  pairs  found  in  the  above  process  are 
recorded  as  the  subscripted  variable  XX,   beginning  with  XXfi.     The 
number  of  "loop  pairs"  is  recorded  as  Z,    (M1365). 

To  determine  the  "loops  not  touching  each  other,   taken  3  at  a 
time"  (hereafter  referred  to  as  "triple  products"  or  "loop  triples"), 
it  is  now  easier  to  use  the  loop  numbers  of  the  "loop  pairs".     These 
pairs  are  recorded  as  XXQ  to  XX?„ 

The  procedure  is  more  readily  explained  by  an  example.     Con- 
sider 5  loop  pairs  recorded  as:" 

(1)  (2)  (3)  (4)  (5) 

XXn    y\X1  XX,,    XXn  XX,     XXr  XX„    XX_  XXn  XX. 

01  2  6  45  67  89 

AC  B  B+l  D  D+l 

The  variables  A,   B,    C,   D  are  defined  initially  as  the  XX  subscripts 
as  indicated.     In  order  to  have  any  "triple  products,"  one  number  of 
each  oT  two  pairs  must  be  the  same.     By  the  manner  in  which  the  pairs 
were  recorded,    the  first  number  of  the  pair  is  always  less  than  the 
second  and  is  always  equal  to  or  less  than  the  first  number  of  each 
succeeding  pair.     Therefore,   in  order  for  loop  number  XX.    to  be 
one  of  a  loop  triple,   XXR  must  equal  XX..     If  not,    then  all  four  vari- 
ables are  increased  by  two  and  the  'check  is  made  agin.     If  XX.    =  XXRJ 
in  order  to  have  aUoop  triple,  loops  XX ~  and  XXR   .   must  not  touch 
each  other,   and  hence  will  be  recorded  somewhere  farther  on  as  a 
pair.     Check  for  this  by  using  XX^..     For  this  to  occur  XXn  =  XX^,    if 
not  then  D  is  increased  by  two  and  checked  again  (until  D  =  2Z, ).     If 
XX_j  =  XXp,    then  to  have  a  loop  triple  XX~    ..   must  equal  XXR   .  .     If 
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this  is  found  to  be  true,    the  loop  numbers  of  the  three  non-touching 
loops  are  XX.,   XX _  and  XXR   .  .     These  three  numbers  are  record- 
ed as  subscripted  XX  variables  beginning  with  XX?7    .     The  number 
of  triple  products  is  recorded  as  Z2  (M1570). 

The  "loops  not  touching  each  other,    taken  four  at  a  time" 
(hereafter  referred  to  as  "4' s"  or  "products  of  four")  are  found  by 
an  extension  of  the  above  process. 

Consider  the  same  5  loops  pairs  as  above,   and  two  triple  pro- 
ducts: 

(1)  (2)  (3)  (4)  (5) 

XXn  XX.  X-Xq  XXn  XX*  AAr  AAn  2\2\ri  Q  Q 

A         C  B  D 


(1)  (2) 

XX10    XX11    XX12  XX13    XX14    XX15 


For  a  "product  of  four"  to  exist  XXA  =  XXg  =  XXD  and  XXCXXB+1 
XXn   -     do  not  touch  each  other.     If  the  latter  statement  is  true,    then 
XX~,   XX,,.,,   XXni1   have  been  recorded  as  a  triple  product  and  the 
program  follows  a  procedure  similar  to  that  described  for  the  triple 
products.     If  any  "4' s"  are  found,    the  loop  numbers  are  XX.,   XX~, 
XXR    ..   XX  and  these  are  recorded  as  XX  variables  immediately 

following  the  loop  triples.     The  first  loop  number  of  the  "4' s"  would 
be  XXQ„      ,  Q7    .     The  number  of  "4' s"  is  recorded  as  Z,  (M1795). 

3.7.     Loops  and  Products  of  Loops  not  Touching  Paths  (M1805-2500) 
The  same  reasoning  regarding  the  extent  of  loop  products  rou- 
tines applies  equally  to  this  subroutine  and  the  physical  limitation  is 
placed  at  having  no  more  than  products  of  three  loops  not  touching  one 
another  and  not  touching  a  path.     This  is  not  considered  a  limiting  fac- 
tor since  the  complexity  of  the  system  with  greater  than  this  number 
of  loop  products  not  touching  paths  would  be  such  that  other  limitations 
previously  mentioned  would  prevail. 
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The  method  used  here  is  similar  to  that  of  the  L     subroutine. 

r 

The  procedure  is  to  find  all  loops  not  touching  path  1,    then  all  loop 
pairs  and  loop  triples  not  touching  path  1.     After  these  are  stored, 
the  process  starts  again  with  path  2  and  so  on  until  completion. 
Consider  the  following  as  having  been  previously  recorded: 

2  direct  paths; 


(1) 
YY 


1 


YY< 


(2) 
YY 


A+2 


YYA+1     YYB+f° 


(1) 


4  loops 
(2)  (3) 


(4) 


ZZC+1      ZZD+1      ZZE+1 


zz„=o  ZZ    =0    ZZ    =0     ZZ    =0 
C  C  E  F 


3  pairs 


loops 

1,2 

1,4 

2,4 

1  triple 

loops 
1,2,4 


no  "4«b" 


The  pairs  and  triples  are  stored,   as  noted  before,   as 


XX«      XX.       XX„       XX„       XX  i       XXp       XXp      XXn      XXq 


and  the  following  variables  have  been  previously  defined: 
Y0        number  of  paths  =    2 

Z0        number  of  loops  =    4 

Z|         number  of  pairs  =    3 

Z?        number  of  triples  =     1 

Z0        number  of  "4' s"  =    0 
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G,  first  ZZ  number  of  loop  1  =1 

G2  first  ZZ  number  of  loop  2  =   C+l 

Gg  first  ZZ  number  of  loop  3  =    D+l 

G4  first  ZZ  number  of  loop  4  =    E+l 

The  procedure  is  to  check  the  first  number  of  path  1  against  all  numbers 

of  loop  1,    then  the  second  number  of  path  1,   etc.,    until  YY»+1   (0)  is 

found.     If  a  number  is  found  in  the  loop  equal  to  the  check  number  of 

the  path,    the  loop  touches  the  path  (Ml 855),    the  loop  number  is  in- 

i 
creased  by  1  and  the  check  number  of  the  path  reset  to  YY..     If  all 

mumbers  of  the  path  are  checked  against  a  loop  and  no  common  point 
found,    the  loop  does  not  touch  the  path  and  the  loop  number  is  record- 
ed as  the  next  subscripted  XX,   in  this  example  the  first  loop  found  not 
to  touch  path  1  is  recorded  as  XXg  (Ml 900).     After  path  1  has  been 
checked  against  all  loops  (M1835),   the  number  of  loops  not  touching 
path  1  is  recorded  as  Y,    (Ml  925).     In  the  example,    say  that  loop  1,    3 
and  4  are  found  not  to  touch  path  1.     These  are  recorded  as  XXq=l, 
XX1Q=3,   XXU  =  4,   Yt  =  3. 

To  find  the  pairs  of  loops  not  touching  path  1  and  not  touching 
each  other,   it  is  necessary  only  to  check  the  loop  numbers  of  single 
loops  not  touching  path  1,   against  the  loop  pairs,    recorded  as  XX0  to 

XX 

2Z^.     This  procedure  is  indential  to  that  used  in  the  loop  products 

routine: 

(1)         (2)  (1)        (4)  (2)        (4)  (1)       (2)         (4) 

XXq    XX-        XXp    XXo       -^-^4    XX,-        XXp    XX„    XXp 


(1)         (3)  (4) 

XXg     XX1Q     XXn      XX12  (not  yet  recorded) 

Q  C  B  E 

The  variables  A,   B,   E,   Q  are  defined  initially  as  the  XX  subscripts 
indicated  and  the  recorded  values  of  the  XX  variables  are  as  indicated 
in  parenthesis.     The  procedure  is  to  check  XXQ  against  XX.,  if    not 
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equal,   increase  A  by  2  and  check  again.     In  the  example  XXQ  =  XX.  . 
The  variable  C  then  is  defined  as  Q+l  (=10)  and  XX  ^  is  checked  against 
XX.    , .     These  are  not  equal  hence  loops  (1)  and  (2)  are  not  a  "pair 
not  touching  path  1".     By  the  previously  described  numerical  order  of 
recording  the  XX  variables,    it  is  seen  that  if  XX~  is  greater  than 
XX.    .,   then  further  checks  against  XX.   are  necessary  (M2035).     The 
next  pair  is  checked  in  the  same  manner.     Again  XX    =  XX.    (A  now 
equals  2)  and  XX~  is  checked  against  XX.        (M2025)  and  found  unequal. 
Now,   however,   XX  „,  is  less  than  XX  .        so  the  possibility  remains 
that  XX?  and  XX~  are  a  pair  not  touching  this  path.     C  is  then  increas- 
ed by  one  (M2045)  (to  11)  and  XX^  is  again  checked  against  XX.    ... 
This  time  the  two  are  found  to  be  equal  and  XX~  and  XX~  are  record- 
ed as  XX.  2  and  XX-,.     Now  A  is  increased  by  2  (to  4)  and  XXQ  is 
checked  against  XX.  (Ml 965),   found  to  be  unequal  and  furthermore, 
XX.  is  greater  than  XX„,   so'no  further  test  is  necessary  using  XX„. 
The  procedure  continues  and  no  more  pairs  are  found.     The  number  of 
pairs  of  loops  not  touching  path  1  as  recorded  as  Y2  (M2100).     Since 
only  one  pair  has  been  recorded  there  can  be  no  "triples  not  touching 
path  l".     If  more  than  one  pair  had  been  found,    these  would  have  been 
checked  for  triples  in  a  similar  manner  to  that  used  for  loops  and  the 
number  of  triples  is  recorded  as  Y,  (M2335).     If  none  are  found,    Y^ 
is  recorded  as  zero. 

The  next  step  is  to  determine  the  end  of  the  present  path  in  the 
YY  list  (M2385,    2390),    check  to  see  if  there  are  any  more  paths  (M2400) 
and  if  so,    repeat  the  process  for  the  next  path. 

The  entire  XX  list  is  summarized  in  the  next  paragraph.     In  ad- 
dition to  the  XX  variables,    the  following  have  been  defined  by  this  sub- 
routine: 

Y1         Number  of  loops  not  touching  path  1 

Y?        Number  of  pairs  of  loops  not  touching  path  1 

Yo        Number  of  loop  triples  not  touching  path  1 

Y4 

5    {   Same  as  Y. ,    Y?,    Y,,  applied  to  path  2 
Y6 
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Y3Y0-2 

3Y«-1     {Same  as  above  applied  to  the  last  path 

Y 
3Y0 


3.8.     Performance  Function  Subroutine  (M2500-3560) 

Before  discussing  this  subroutine,   it  is  best  to  summarize  the 
XX  variables  recorded  in  the  routines  of  paragraphs  3.  6. ,   and  3.  7: 


XXQ -      Loop  Pairs XX2, 


prl 


XX2Z -Loop  Triples-  -  -   XX2Z+3Z2-1 

TK2zi¥3:z^  ------  Lo°p  "4's" "aZj+aZj-rtZj-i 

Define  E  =  2ZX+3Z2+4Z3 

XXp,  —  —    Loops  not  touching  P.  -  -  XX^.y      i 

XXE+Y   ~  ~  "LooP  Pairs  not  touching  Pl  -  -   XX£+Y  +2y   _1 

XXE+Y1+2Y2  TriPleS  not  touching  Pl XXE+Y1+2Y2+3Y3-1 

Define  F  =  E  +  Y1+2Y2+3Yg 


XX 


-  -  -  -   Loops  not  touching  P2 


etc.  for  all  paths. 
The  last  number  recorded  is  XX„.   where  G  is  defined  as: 

G  =  2Z1+3Z2+4Z3+Y1+2Y2+3Y3+Y4  +  -  -  -  +3Y3y 

As  stated  in  Chapter  2. ,    this  routine  is  written  to  provide  a  per- 
formance function  of  the  form 

G1G2G3  +  G4G5  + 


1  +  G  G,  G    +  G  ,G    +  - 

a    b    c         d    e 
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The  results  of  the  routine  are  the  subscripts  of  the  G'  s  with  appro- 
priate algebraic  signs.     This  routine  is  by-passed  if  not  desired 
(M2501).     The  by-pass  is  effected  as  stated  in  paragraph  3.2.,   by 
entering  a  zero  for  the  variable  SS  on  the  first  data  card. 

The  subroutine  is  quite  naturally  divided  into  two  parts,    namely, 
a  Numerator  and  Denominator.     Recall  that  the  numerator  is: 

k    L  rk  rk  rk  J 

In  order  to  punch  this  function,   it  is  necessary  only  to  know 
where  the  particular  values  are  stored  in  the  computer.     The  method 
is  very  straight  forward  and  will  be  described  in  detail  only  where 
particular  sections  of  the  program  might  require  clarification. 

The  YY  and  ZZ  lists,    of  path  and  loop  components  respectively, 
also  contain  the  summers,   which  are  not  desired  in  the  equation. 
Hence,  any  number  in  these  lists  greater  thaVi  the  number  of  compon- 
ents is  not  punched  (M2600). 

For  convenience,   the  first  YY  number  of  each  path  is  recorded 
(M2515  and  M2635)  in  the  same  manner  as  the  first  ZZ  number  of 
each  loop  was  recorded  earlier.     All  components  (excluding  summers) 
of  path  1  are  punched  successively  (M2600,    2605)  with  no  sign  of  posi- 
tive and  with  a  minus  sign  of  the  term  of  negative  (M2585,   2595).     The 
second  term  of  the  numerator  (P,    ZL  J  s  P-  L-.  +P.  L?1   ...   is 
a  series  of  Y.   terms,    the  first  of  which  contains  all  the  components 
of  path  1  and  all  the  components  of  the  first  loop  not  touching  path  1, 
This  loop  number  is  found  in  the  XX  list,    the  first  ZZ  number  of  the 
loop  is  known  so  the  components  of  the  loop  are  located  and  the  entire 
term  is  punched. 

Throughout  the  routine,   there  are  many  occassions  to  use  a  par- 
ticular group  of  equations   or  commands,    equations  locating  the  number 
in  the  XX  list,   punch  plus  and  punch  minus  commands,   etc.     To  obviate 
the  necessity  of  Writing  these  same  equations  many  times,    two  "routine 
indicators"  P  and  R,   are  used  to  transfer  control  back  to  the  desired 
program  locations. 
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Continuing  the  numerator,    the  next  series  of  terms  is 


Pk    2   j|  Lrk 


and  the  procedure  is  similar  to  that  used  before.     After  all  terms 
have  been  punched,   the  same  process  is  carried  out  for  path  2  etc. 
until  the  entire  numerator  has  been  punched.     To  fully  understand  the 
mechanics  of  the  program  requires  a  step  by  step  analysis  of  a  sam- 
ple system.     This  is  done  for  a  sample  system  of  Chapter  4. 

After  all  numerator  terms  have  been  punched,  the  program 
shifts  to  the  denominator  (M3083).  Hence,  the  routine  is  simpler 
since  the  denominator  of  the  performance  function  equation  is 

1  +  2L     +    2  II  L.      +  2  II  L.      +  2  II  L     +.. 
,r  2      r  3       r  4       r 

The  procedure,  is  the  same  as  for  the  loops  of  the  numerator  and  the 
components  of  each  loop,   loop  pair,   etc.   are  punched  successively 
with  appropriate  algebraic  signs. 

3.  9.     Frequency  Response  Calculation  (M3570-M5285) 

Having  set  up  all  terms  of  the  performance  function,   there  re- 
mains only  the  task  of  devising  a  scheme  for  entering  the  transfer 
functions  of  the  system  components,    substitute  S  =  j  u,   assign  a  fre- 
quency value  and  combine  the  results  into  the  terms  of  the  perform- 
ance function.     As  described  in  Chapter  2,    the  comppnent  transfer 

functions  are  a  ratio,    the  numerator  and  denominator  of  which  are 

2 
comprised  of  products  of  terms  of  the  form  AS     +  BS  +  C,   of  which 

any  of  the  coefficients  may  be  zero.     It  is  desired  to  enter  on  data  cards 
only  the  coefficients.     Since  each  transfer  function  may  contain  any 
number  of  terms,   provisions  must  be  made  for  recognizing  the  number 
of  terms  for  each  component.     This  is  accomplished  by  means  of  an 
array  of  the  subscripted  variable  G.     For  ease  of  entering  the  infor- 
mation a  complete  form  is  included  with  the  sample  problem  of  Chapter 
4. 
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To  understand  the  form  in  which  the  data  is  entered  and  the  re- 
quired information  obtained,   a  simple  example  is  shown  in  Fig.    3.2. 

As  stated  previously,   in  order  for  the  program  to  be  properly 
executed,   a  scheme  must  be  established  whereby  the  terms  of  the 
array  are  recognized.     The  scheme  chosen  to  accomplish  this  is  to 
enter  the  G  subscript  of  the  first  number  of  the  last  term  of  each 
numerator  and  denominator.     These  numbers  will  be  referred  to  as 
"numerator  stops"  and  "denominator  stops"  respectively.     The  num- 
bers are  shown  circled  in  Fig.   3.2.   and  are  entered  in  the  following 
order: 

(1)  Numerator  stops  listed  consecutively,    6  to  a  card 

(2)  Denominator  stops  listed  consecutively,    6  per  card. 
The  last  card  of  each  may  not  necessarily  be  full.     In  the  example  of 
Fig.   3.2.  ,   only  one  card  is  needed  for  each.     The  first  card  would 
contain:    0,    15,    24  and  the  second  card:    6,   21,   27.     These  numbers 
are  entered  into  the  G- array  immediately  following  the  last  number 
of  the  transfer  function  array,   with  a  blank  space  between  the  two 
sets  of  numbers.     This  blank  space  is  given  the  value  of  -3,   which 
is  necessary  for  manipulations  later  on. 

Example  of  the  Transfer  Function  Array 

Consider  a  system  of  3  components  with  transfer  functions: 


1) 


2) 


3) 


S  +  1 


S    (5S  +  1) 

62  (.2S  +  1)  (3S  4   1) 
(4S  +  1)  (S2  +  7S  +  1) 


2S  +  1 


The  coefficients  are  entered  in  the  following  order: 

1.  All  numerator  terms  of  component  1 

2.  All  denominator  terms  of  component  1 

3.  All  numerator  terms  of  component  2 

4.  All  denominator  terms  of  component  2 

5.  All  numerator  terms  of  component  3 

6.  All  denominator  terms  of  component  3 

-39- 


The  numbers  are  entered  consecutively,    6  per  data  card.     These  be- 
come a  G-array  as  (numbers  in  parenthesis  and  circled  are  sub- 
scripts of  the  variable  G): 


Numerator  1 


Denominator  1 


Numerator  2 


Constant 
0 


Denominator  2       (18)  1 

I 


21 


Numerator  3  (  24 


Denominator  3      (  27)  1 


1 


Coeff .   of  S 
(1)      1 


(19)      4 
(22)      7 


(25)       0 


(28)       2 


Coeff.    of  S' 
(2)       0 


(3)          0 

to  > 

(4)      1 
(7)      5 

(5) 
(8) 

0 
0 

(9)       62 
(12)         1 

0       x 

(10)      0 
(13)    .2 
(16)      3 

(11) 
(14) 
(17) 

0 
0 
0 

(20)        0 
(23)        1 


(26)        0 


29)        0 


Fig.    3.2.     G- Array  of  Coefficients  of  Transfer  Functions 


One  other  card  is  required  to  be  entered  and  it  contains  the  fre- 
quency information  as:    (M3570) 

Initial  frequency  (radians  /second)   -----  VV 

Frequency  increment  desired    -------     DW 

Final  frequency  (radians/ second)   -  -  —  -     FW 
Since  a  MAC  READ  order  for  reading  in  arrays  of  numbers  must  con- 
tain some  information  as  to  the  number  of  cards  to  be  read,    it  is  neces- 
sary that  the  last  G  subscript  of  the  transfer  function  array  be  entered. 
This  number  (shown  circled  in  Fig.    3.2.    as  29)  is  entered  on  the  fre- 
quency data  card  and  recognized  by  the  computer  as  RR.    (M3570). 
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Having  read  in  the  necessary  data,    the  program  picks  up  the 
transfer  function  of  the  first  component.     Letting  S  =  j  u  and  sub- 
stituting the  initial  frequency  (VV),   a  magnitude  and  angle  is  com- 
puted for  the  component.     Then  the  same  process  is  carred  through 
for  each  component.     Referring  to  Appendix  B,    card  M3640,    the  pro- 
cedure is  as  follows: 

1.  Define  the  variable  H  initially  as  the  subscript  of  the 
G  variable  left  blank  between  the  numerator  and  de- 
nominator stops. 

2.  Establish  the  variables  K  and  L  as  the  limits  in  the 
G- array  of  all  terms  (numerator  and  denominator) 
of  component  1 . 

Space  available  as  subscripted  XX  is  used  as  temporary  storage 
for  the  computation  of  the  magnitude  and  angle  of  the  components.     The 
magnitude  and  angle  of  each  term  of  component  1  is  computed  by  M3690 
and  M3730.     The  command  DO  TO  is  used  to  accomplish  this  and  means 
to  carry  out  the  operation  specified  for  various  values  of  the  variable 
K  from  its  present  value  in  steps  of  (3)  to  a  final  value  L.     The  opera- 
tion is  merely  to  determine  the  square  root  of    the  sum  of  the  squares 
of  the  complex  number  representing  each  term  of  component  number  1. 
The  angle  is  found  by  M3730  as  the  ARCSIN  of  the  imaginary  part  of 
the  complex  number  divided  by  the  magnitude.     The  ARCSIN  function 
will  give  results  between!   /  2,   hence,    it  is  necessary  to  make  the  check 
of  M37  60  which  states  that  if  the  real  part  of  the  complex  number  is 
negative,   subtract  the  angle  from  n.     If  the  real  part  is  positive,    the 
operation  is  ignored. 

To  see  the  manipulations  more  clearly,    refer  to  Fig.    3.2. 

R  =  29 

H  =  33 

K  =  G33  +  3  =  0 

L  =  G34  =  6 

This  is  seen  to  indicate  3  terms  of  component  1,   and  by  M3690  and 
M3730,   a  magnitude  and  angle  is  computed  for  each  of  them.     To  ob- 
tain the  magnitude  and  angle  of  the  entire  component,    the  numerator  and 
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denominator  terms  are  separated,   each  manipulated  as  necessary, 
then  the  magnitudes  divided  and  the  angles  subtracted  to  give  the 
results.     This  is  accomplished  by  M3780  to  M3925.     Referring  again 
to  Fig.   3.2.,   and  starting  with  M3 7 85: 
K  =  0 

M  =  G29  +  1  =  ° 
This  means  that  there  is  only  one  term  in  the  numerator  and  the 
numerator  magnitude  then  is  the  magnitude  already  computed  for  that 
term.     The  numerator  angle  is  found  by  adding  the  angles  of  all 
numerator  terms  (M3835,    M3840)  but  in  this  case,    there  is  only  one. 

A  similar  procedure  is  carried  through  for  the  denominator 
terms : 

(M3855)    K  =  G2g+3  =  3  established  S  as  the  G  subscript 

of  the  first  denominator  term.     The  two  denominator  terms 

are  multiplied  by  M3865  and  3870,   in  two  steps,    K  =  3,  6. 

M3850      C  =  1 

M3870       C  =  (1)    (XX^+3)  =  XXE+3 

M3870       C  =  (XXE+3)    (XXE+6) 

In  a  similar  manner,   the  2  denominator  angles  are  added.    (M3895, 
M3900) 

It  is  desired  now  to  compute  a  magnitude  and  angle  for  the  cora- 
ponent  as  a  subscripted  variable  Gf  in  the  interests  of  space  economy, 
since  later  subroutines  use  the  XX  variable.     Therefore,    F  is  com- 
puted as  the  next  available  G  number,    and  the  component  magnitude 
is  computed  by  dividing  numerator  by  denominator  and  identified  as 
G-p.     It  is  desired  to  have  all  component  magnitudes  in  successive  G 
numbers  followed  by  all  component  phase  angles.     Hence,   with  Q 
components  the  phase  angle  of  component  one  is  computed  by  subtract- 
ing denominator  angle  from  numerator  angle  and  identifying  the  result! 

as  GF+Q' 

The  procedure  now  is  to  return  to  the  beginning  of  this  subroutine 

and  carry  out  the  above  process  on  the  remaining  components  in  suc- 
cession.    Since  the  component  magnitude  and  phase  angles  will  be  used 
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later,    it  is  convenient  to  define  the  variable  J  (M3970)  as  the  last  G 
number  preceding  the  first  component  magnitude.     Then  to  find  the 
magnitude  of  component  N,   the  location  is  known  to  be  Grj,^. 

It  would  be  satisfactory  now  to  substitute  these  components 
magnitudes  and  phase  angles  into  the  performance  function  as  punch- 
ed in  paragraph  3.  8.  ,   perform  the  required  algebra  and  obtain  the 
response  of  the  system  for  this  frequency.     However,   the  perform- 
ance function  is  not  retained  as  such  and  calculation  of  the  perform- 
ance function  is  desired  to  be  optional.     Since  the  XX  variables  are 
stored  as  loop  numbers,   each  loop  number  appearing  many  times,   a 
more  convenient  and  faster  method  is  to  compute  the  magnitude  and 
phase  angle  of  each  loop.     Then  the  loop  response  is  used  in  comput- 
ing the  response  of  the  loop  product  terms.     All  of  these  are  terms 
in  the  performance  function  and  will  be  referred  to  herein  as  "term 
responses". 

The  first  step  is  to  devise  a  storage  scheme  for  these  term 
responses  so  they  may  be  easily  combined.     The  denominator  of  the 
performance  function  is  computed  first  and  the  storage  scheme  used 
is  (M4005  -  4025):     First  available  XX  variable:     XX£ 

XXE+1"  "  XXK  XXK+1  XXL 

Loop  magnitudes  Loop  phase  angles 


xxL+r  -  -  -  -  -xxc       XXC+1 

Real  part  of  all  Imaginary  part  of 

term  responses  all  term  responses. 

The  denominator  is  a  summation     of  term  responses: 

(1  +    2  Lr  +  2  n  Lr  +  .    .    „) 
2 

After  computing  these  term  responses  and  converting  to  rectangular 

form,   all  that  is  necessary  is  to  add  the  real  parts  (XX.     ,   to  XX„) 

and  the  imaginary  parts  (XX„ , .   to  XX^,,,,         „         „         ,-,    \     ,     _ „«„* 

6  J  r  C+l  C+Zn  +  Z1   +  Zp  +  ZJ,    to  convert 
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the  result   to  polar  form  and  retain  it  until  after  the  numerator  res- 
ponse is  computed. 

In  computing  loop  magnitudes  recall  that  the  ZZ  list  of  loop  com- 
ponents (paragraph  3.4.)  has  a  sign  determined  by  the  number  of  minus 
signs  encountered.     For  example,   a  single  loop  is  recorded  as: 


+ 


zzx  = 

-3 

zz2- 

1 

zz3  = 

2 

zz4  = 

0 

3 

7) 

■ 

JL 

F 

2 

As  explained  in  Chapter  2  and  from  elementary  feedback  theory 
the  sign  of  this  loop  is  (+).     Hence,    an  initial  multiplier  of  -1  is  used 
(M4045)  and  an  initial  phase  angle  of  0  (M4050).     Component  magnitudes 
and  phase  angles  have  been  Stored  as  Gy,  .   and  Gj+a+q  respectively. 
The  loop  magnitude  and  phase  angle  is  computed  (M4055  -  4115)  using 
only  components  (summers  are  excluded),   with  the  appropriate  sign, 
and  recorded  as  described  in  the  storage  scheme  above.     It  is  required 
that  the  loop  response  be  retained  in  both  polar  and  rectangular  form 
(M4130  -  4160). 

The  response  of  the  terms  representing  loop  pairs  is  obtained 
by  (M4210  -  4290)  using  the  loop  numbers  of  the  pairs  from  XX  storage 
beginning  with  XXQ  and  multiplying  magnitudes  and  adding  phase  angles 
of  the  loops.     Since  these  loop  pair  responses  are  never  multiplied  in 
the  performance  function  equation,    they  are  converted  to  rectangular 
form  (M4250  -  42  60)  and  retained  only  in  this  form. 

By  exactly  the  same  method,    the  term  responses  for  the  loop 
"triples"  and  "4'  s"  are  computed  and  retained  in  rectangular  form 
(M4310  -  4500).     The  denominator  response  is  then  obtained  by  adding 
all  real  parts  of  the  term  responses,    stored  consecutively  as  XX,     « 
to  XX„,   and  all  the  imaginary  parts,   XX_       to  XX„  .     T  is  a 

variable  used  to  increase  the  XX  variable  by  one  each  time  a  new  term 
was  added,   hence  it  can  be  used  as  a  check  to  determine  the  end  of  the 
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list  to  be  added.     There  are  T-l  numbers  in  each  list.     The  result  of 
the  additions,   a  complex  number,   is  converted  to  polar  form  (M4580- 
4630)  and  retained  as  denominator  magnitude,    (DM)  and  denominator 
phase  angle  (DA). 

Note,    that  the  conversion  from  rectangular  to  polar  form  is  a 
routine  that  is  to  be  used  more  than  once,   so  the  routing  identifier  (X) 
is  used  (M4570).     As  noted  previously,    the  ARCSIN  function  given  only 
angles  between  +     /  2,    so  if  the  summation  of  the  real  part  is  negative 
the  correct  angle  must  be  computed  (M4605). 

Obtaining  the  numerator  response  is  only  slightly  more  com- 
plicated since  it  is  a  summation  of  terms  multiplied  by  a  path  function 
and  then  the  products  are  summed: 

ZPk    t    1+SLrk+2^    Lrk+---    ] 

The  procedure  is  to  determine  the  path  response  (Pk)  for  path  1,    in 
polar  form,    then  determine  the  term  response  of  each  of  the  terms 
indicated  in  the  equation  that  are  associated  with  this  path,   perform 
the  addition  indicated  in  the  brackets  and  convert  to  polar  form.     This 
result  is  then  multiplied  by  the  path  response,   the  final  result  convert- 
ed to  rectangular  form  and  stored. 

The  procedure  is  then  repeated  for  the  second  path,   and  so  on. 
After  all  paths  have  been  completed  the  i  ihal  summation,  Z  P,    [  1  +  .  .  ] 
is  performed  and  the  result  converted  to  polar  form. 

The  loop  responses  previously  computed  are  needed  here  but  the 
other  denominator  terms  are  not,   hence,    the  same  storage  may  be  used 
for  the  numerator  terms.     In  the  program  of  Appendix  B,   this  is  only 

partially  true,   in  that  it  was  chosen  to  use  only  the  space  after  XXn   7    : 

^0 
(B  =  C  +  ZQ) 

XXE+1  "  "       XXL+1  XXC+1 

Loop  magnitude  Real  part  of  all  Imaginary  part 

and  phase  angle  denominator  term       of  loop  responses 

responses 


XXB+1  ------  xxH+l 

Real  part  of  all  numerator  Imaginary  part  of  all  numerator 

term  responses.  term  responses. 
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The  magnitude  and  phase  angle  of  the  fir^t  path  are  computed 
(M4690-  M4760)  in  the  same  manner  as  that  of  the  loops  except  that 
the  sign  of  a  path  function  is  minus  if  the  path  in  the  YY  list  has  a 
minus  sign.     The  path  magnitude  and  phase  angle  are  retained  tempor- 
arily as  V  and  W  respectively.     (M4720  -  4730). 

The  numerator  terms  associated  with  the  path  are  located  in  XX 
storage  from  the  stbrage  scheme  summarized  in  paragraph  3.8.     The 
first  group  of  these  terms  (2  L   .  )  is  a  summation  of  loops  for  which 
the  real  and  imaginary  parts  of  their  response  have  been  previously 
computed,    so  need  only  be  rewritten  as  the  appropriate  XX  variable 
(M4770  -  4825). 

The  second  group  of  numerator  terms  (ZHL   ,  )  and  the  third 
group  (S  n    L   .  )  are  determined  exactly  as  the  S  II  L     and  ZTn  L 
terms  of  the  denominator  (M4835  -  4915  and  M4925  -  5015). 

Having  determined  the -response  of  all  terms  associated  with 
this  path,   their  real  parts  +  1  and  imaginary  parts  are  summed  (M5030- 
5085).     The  routing  identifier  X  =  1  is  used  and  the  program  jumps  to 
the  routine  for  converting  rectangular  to  polar  form  (M5095). 

Up  to  this  point,    the  following  have  been  obtained:    V/_W_  [  P/  A 
(for  1  +  2  L   ,    +  .    .    )].     The  response  for  all  terms  associated  with 
this  path  is  then  computed  (M5100  -  5120)  and  retained  as  an  XX 
variable.     If  any  more  paths  exist  (M5130)  the  same  routine  is  carried 
out  for  each. 

After  all  paths  have  been  completed,    the  real  and  imaginary  parts 
are  summed  (M5155  -  5195).     The  storage  locations  of  these  values  are: 

XXC  -  -  -   I  XXB<+1  XXB+2   I  XXB+3  XXB+4   I  etc. 
first  path  second  path 

Each  time  a  path  was  completed  B  was  increased  by  2  in  preparation 
for  working  space  of  the  next  path.  After  all  paths  are  completed,  B 
is  set  back  to  its  original  value  (M5156)  to  permit  summing. 
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The  program  now  jumps  to  the  conversion  of  rectangular  to  polar 
form  routine  (M5210),    the  result  of  which  gives  the  numerator  magni- 
tude and  numerator  phase  angle.     The  system  magnitude  ratijp  (M)  and 
phase  angle  (P)  are  found  (M5215  -  5220).     The  final  result  is  desired 
with  the  phase  angle  given  as  a  positive  angle  between  0  and  360  de- 
grees,  hence,    the  two  equations  (M5235  and  5255).     The  possibility 
exists  that  M  might  be  zero  or  the  DM  might  be  zero  (M  would  be  in- 
finite).    In  either  case,   special  punch  orders  are  provided  (M5212,   5275 
and  5237,    5253). 

The  following  information  is  obtained  from  the  punched  result: 
M         Magnitude  Ratio 
P  Phase  Angle  in  Degrees  (0-360) 

AA  Log  M  (MAC  computes  logs  to  the  base  e  so  it  is 

converted  to  the  base  10) 
20AA  Magnitude,  Ratio  in  d.b. 

VV  Frequency  (radians/ sec) 

VV/2I1  Frequency  (cps) 

3.10.     Special  Features 

There  are  certain  special  Situations  that  are  not  strictly  a  part 
of  any  of  the  general  subroutines  of  the  preceding  paragraphs,   but 
might  require  explanation.     One  such  situation  might  arise  when  the 
system  performance  function  is  known  in  factored  form,   perhaps  from 
a  pole-zero  configuration.     In  this  case,    the  block  diagram  of  the  sys- 
tem is: 


Component  1  may  consist  of  any  number  of  terms.     There  are  no  loops 
and  only  one  path.     The  calculation  of  frequency  response  of  the  system 
is  made  treating  the  system  as  one  component  and  the  entire  program 
up  to  and  including  the  performance  function  routine  is  by-passed.     This 
by-pass  is  effected  by  the  presence  of  a  zero  on  the  first  data  card  in 
the  position  designated  as  "number  of  summers"  (M0035,    0036).     The 
number  of  loops  (Zfl)  is  recorded  as  zero  in  the  by-pass  (M2565)  and 
the  variable  E,   used  to  establish  storage  locations  in  the  XX  list  is 
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assigned  a  value  of  zero. 

The  "input  data  array"  is  not  entered.     The  first  card  must  con- 
tain six  numbers  since  the  READ  Order  (M0035)  says  read  six  values. 
However,    the  only  ones  with  any  meaning  in  this  case  are  the  third, 
where  a  zero  is  entered  for  S,    and  the  second  where  the  number  of 
components  is  entered  as  1.     The  remainder  of  the  data  cards  are 
entered  (M3570,    3585,    3;595,    3605)  exactly  as  explained  in  paragraph 
3.  9. 

A  second  special  situation  arises  when,   having  computed  a  fre- 
quency response,   it  is  desired  to  change  the  frequency  range,    fre- 
quency increment,    or  a  transfer  function  of  one  of  the  existing  com- 
ponents.    In  this  case,    it  is  desired  that  all  of  the  loops,   paths,   and 
combinations  leading  to  a  performance  function  be  retained,    and  only 
that  part  of  the  program  pertaining  to  the  calculation  of  frequency  res- 
ponse be  repeated.     A  simple  procedure  for  handling  this  is  to  let  the 
program  run  to  completion,   have  new  data  cards  available,    and  start 
again  at  M3570.     The  program  is  set  up  in  such  a  manner  that  all 
that  is  required  is  to  push  the  start  switch  and  the  program  will  start 
at  M3570. 

If  in  addition  to  the  above,   it  is  also  desired  to  insert  a  new 
component  or  change  the  configuration  of  the  existing  circuit,   it  is 
still  not  necessary  to  read  the  entire  load  deck  again.     A  "program 
start"  may  be  entered  directly  into  the  computer  in  basic  machine 
language.     This  entails  setting  up  on  an  "index  register",   the  address 
of  the  appropriate  command  in  the  computer.     This  is  a  minor  task 
but  to  ensure  proper  operation,   a  detailed  step-by-step  procedure 
follows  : 


Detailed  Step-by-Step  "Re-run"  Procedure 

1.  Prepare  data  cards  as  in  original  situation. 

2.  Set  1698  in  the  ADDRESS  SELECTION  SWITCHES  on  the  main 
console.     These  switches  are  plainly  labeled. 

3.  Set  computer  CONTROL  SWITCH  TO  MANUAL. 
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4.  Press  PROGRAM  RESET  key. 

5.  Press  TRANSFER  key. 

6.  Set  CONTROL  SWITCH  to  RUN, 

7.  Press  PROGRAM  START  key. 
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CHAPTER  4 
SUMMARY  OF  PROCEDURES  FOR  USING  THE  PROGRAM 


4.  1.     General 

In  Chapter  2,    the  general  plan  of  the  program  was  discussed 
with  the  intent  of  describing  its  application  to  any  computer.     Chapter 
3  covered  the  details  of  the  program,   of  Appendix  B,   written  in  MAC 
language  for  the  IBM  650  and~its  associated  equipment.     The  mechanics 
of  the  program  and  the  method  used  to  obtain  the  frequency  response 
have  been  adequately  discussed.     In  this  chapter,   a  summary  of  the 
procedures  required  by  the  user  will  be  outlined  with  application  to  a 
specific  example. 

It  is  noted  that  MAC  uses  plain  language  for  the  basic  program 
but  that  this  program  is  then  assembled  in  "machine  language"  prior 
to  operation.     The  "load  deck"  referred  to  in  Chapter  2  is  a  machine 
language  version  of  the  MAC  program.     Any  other  information  entered 
in  the  form  of  data  cards  must  be  in  a  form  recognized  by  the  computer 
and  independent  of  MAC.     Data  cards  contain  only  10-digit  numbers, 
six  or  less  to  a  card,    and  the  numbers  are  in  "floating  point"  form. 
Floating  point  is  merely  a  method  of  placing  the  decimal  point  in  a 
number  always  written  with  the  first  digit  as  non-zero.     The  first  8 
digits  are  the  numerical  value  and  the  last  two  locate  the  decimal  point. 
For  a  number  with  the  decimal  point  at  the  extreme  left  (e.  g.     .  62), 
the  9th  and  10th  digits  are  50.     For  each  digit  that  the  decimal  point 
is  located  to  the  right  or  left  of  the  first  non-zero  digit  of  a  number, 
1  is  added  or  subtracted  from  50.     For  example: 


62.0 

is  written  in  floating  point  as 

62000000  52 

6.20 

as 

62000000  51 

.620 

as 

62000000  50 

.062 

as 
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62000000  49 

4.2.     Data  Cards  Required   ' 

The  data  cards  are  placed  in  the  card  reader  feed  directly  on 
top  of  the  MAC  or  load  deck  and  are  read  at  the  appropriate  time  by 
a  READ  order  in  the  program.     The  number  of  data  cards  required 
will  vary  with  the  complexity  of  the  physical  system  being  analyzed. 
However,    the  type  of  data  falls  into  six  general  groups.     A  new  card 
must  be  punched  at  the  beginning  of  each  of  these  groups  since  they 
are  read  in  by  different  READ  orders.     The  last  card  of  each  group 
is  not  necessarily  filled.     The  six  general  groups  and  the  data  cards 
of  each  are  as  follows: 

(a).  General  Circuit  Information    (1  card) 

(1)  Maximum  number  of  inputs  to  any  unit K 

(2)  Number  of  component  boxes  in  the  block 
diagram  of  the  system ^Q 

(3)  Number  of  summers S 

(4)  Number  of  the  input  unit CC 

(5)  Number  of  the  output, unit DD 

(6)  Is  Performance  Function  desired?    ...      SS  =  0  if  No 

1  if  Yes 

(b)  Input  Data  Array  (M  cards) 

(where  M  =  QQ  +  S  +  1,   unless  K  is  greater  than  4,    in 
which  case,    the  number  of  values  on  each  card  exceeds 
6  and  additional  cards  must  be  punched) 
The  "input  data  array"  described  in  paragraph  2.2.   and 
3.3.    is  entered,    one  row  per  card  (unless  the  row  con- 
tains more  than  6  numbers,    in  which  case,    an  additional 
card  is  used  for  each  multiple  of  6).     Each  row  of  the 
array  is  started  on  a  new  card.     The  last  row  of  the  array 
and  hence  the  last  card  contains  all  zeros. 

(c)  Frequency  Data    (1  card) 

(1)  The  lower  limit  of  the  frequency  range  over 
which  the  response  is  desired VV 

(2)  The  increment  of  frequency  over  the  fre- 
quency range DW 
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(3)  The    upper  limit  of  the  frequency  range  .    .    .     FW 

(4)  The  last  subscript  number  of  the  G-array 

of  Transfer  Functions RR 

Note  that  only  four  numbers  on  this  card  will  be  read.     The  last 
two  spaces  will  not  and  zeros  may  be  entered  in  these  spaces  if  de- 
sired: 

d.  The  Transfer  Function  Array 

The  number  of  cards  here  will  vary  with  the  complexity 
of  the  physical  system.     Numbers  on  these  cards  are 
punched  in  succession,   six  to  a  card,    until  all  transfer 
function  coefficients  have  been  punched.     As  described 
in  paragraph  3.  10.  ,    the  order  of  the  coefficients  enter- 
ed for  each  term  is:     constant,    coefficient  of  S,   coef- 

2 
ficient  of  S    .     The  order  in  which  the  terms  are  entered 

is: 

All  numerator  terms  of  component  1 

All  denominator  terms  of  component  1 

All  numerator  terms  of  component  2 

All  denominator  terms  of  component  2 

And  so  on  for  each  component;   the  last  card  may 

contain  3  or  6  numbers. 

e.  Num  erator  Stops 

A  new  card  is  started  for  this  group,    regardless  of 
whether  the  last  card  of  group  (d)  is  filled.     Numbers  on 
these  cards,   as  explained  in  paragraph  3.9.,    represents 
the  subscript  of  the  G  number  in  the  transfer  function 
array,    containing  the  first  coefficient  (constant)  of  the 
last  term  of  the  numerator  of  each  component  transfer 
function.     These  numbers  are  punched  consecutively,    6 
to  a  card,   except  that  the  last  card  may  contain  less  than 
six. 
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(f)  Denominator  Stops 

A  new  card  is  started  for  the  first  numbers  of  this  group. 
As  explained  in  paragraph  3.9.,   numbers  on  these  cards 
represent  the  subscript  of  the  G  number  in  the  transfer 
function  array,  containing  the  first  coefficient  (constant) 
of  the  last  term  of  the  denominator  of  each  component 
transfer  function.     These  numbers  are  punched  consecu- 
tively,   6  to  a  card. 

4.3.  Results  of  the  Program 

As  described  in  Chapter  3  and  illustrated  by  the  example  in 
this  chapter,    the  result  of  the  performance  function  routine  (if  such 
is  desired)  is  a  series  of  groups  of  numbers  representing  the  com- 
ponent making  up  a  term  in  the  performance  function  equation.     These 
groups  of  numbers  (terms  of  the  equation)  are  preceded  by  the  appro- 
priate algebraic  sign. 

Results  of  the  frequency  response  calculation  are  given,    for  each 
frequency,   as  a  row  of  six  numbers,   under  proper  heading  and  re- 
presents: 

a.  Magnitude  ratio 

b.  Phase  angle  (in  degrees) 

c.  Log  of  magnitude  ratio 

d.  Magnitude  ratio  in  decibels 

e.  Frequency  (radians/ sec) 

f.  Frequency  (cycles/ sec) 

4.4.  Example 

Consider  the  system  represented  by  the  block  diagram  of  Fig.   4.  1. 
with  transfer  functions  as  indicated.     This  is  a  hypothetical  system 
chosen  of  sufficient  complexity  to  adequately  test  all  parts  of  the  pro- 
gram and  to  illustrate  the  advantage  of  digital  computation  over  graphi- 
cal techniques. 

The  first  step  is  to  list  the  information  required  by  the  single  card 
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Output Inputs 

—J—  |  ■ 

1  8  0                    0                    0 

2  9  0                    0                    0 

3  10  0                    0                    0 

4  11  0                     0                     0 

5  12  0                    0                    0 

6  13  0                     0                     0 

7  8  0                     0                     0 

8  -6000 

9  12  0                     0 

10  2-370 

11  3-4-5  0 

12  4  0                     0                     0 

13  5-680 
0    ►  0                  0                     0                     0 


Fig.   4.2.     Input  Data  Array 
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G.  Subs. 

Const.   ( 

3.  Subs 

Coef.  of  ( 
S 

1.  Subs. 

Coef.  of 

s2 

© 

15 

1 

0 

2 

0 

© 

1 

4 

0 

5 

0 

© 

1 

7 

0 

8 

0 

9 

© 

0 
1 

10 
13 

1 
1 

11 
14 

0 
0 

© 

10 

16 

0 

17 

0 

18 

0 

19 

1 

20 

0 

© 

1 

22 

1 

23 

0 

© 

1 

25 

.2 

26 

0 

© 

1 

28 

1 

29 

1 

© 

5 

31 

0 

32 

0 

© 

1 

34 

0 

35 

0 

© 

1 

37 

0 

38 

0 

© 

1 

40 

.5 

41 

0 

© 

1 

43 

0 

44 

0 

© 

1 

46 

.03 

© 

RR 

0 

Fig.    4.3.     Transfer  Function  Array 
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of  group  (a).     Values  for  the  variables  are  seen  to  be: 

K  =  3 
QQ  =  7 

S  =    6 

CC  =  8 

DD  =  6 

SS  =    1 
i 
Succeeding  steps  in  preparing  data  for  data  cards  are: 

1.  Write  the  input  data  array  as  in  Fig.   4.  2. 

2.  List  the  frequency  data.     For  this  example,    a  frequency 
range  of  .  1  to  6.  0  radians  per  second  was  chosen  with  a 
frequency  increment  of  .4  radians/ sec.     Note,    that  the 
last  item  on  this  card  does  not  pertain  to  frequency  data 
but  is  a  required  entry  and  placed  on  this  card  for  conven- 
ience.    Its  value  will  be  determined  later. 

3.  Write  the  transfer  function  array  as  in  Fig.   4.  3.     Indicate 
in  some  manner  (circles  are  used  in  Fig.   4.  3)  the  numera- 
tor and  denominator  stops  and  the  last  G  subscript  of  the 
array. 

Data  cards  are  now  punched  in  floating  point  form.     These  are  illus- 
trated in  Fig.   4.  4. 

The  program,   followed  by  the  data  cards,   is  then  placed  in  the 
card  reader  and  the  computer  set  in  operation.     Since  the  first  data 
card  indicated  a  performance  function  was  desired,    the  first  result 
is  this  performance  function.     For  this  fairly  complicated  system,   the 
performance  function  is  quite  lengthy  so  is  appended  as  Appendix  C. 

The  frequency  response  results  are  listed  and  plotted  as  shown 
in  Appendix  D.     Note,    that  a  plotter  is  available  and  when  set  to  the 
appropriate  scales  will  produce  a  plot  directly  from  the  punched  cards. 
The  plotter  was  not  used  for  this  example. 

4.5.     Summary  and  Recommendations 

This  report  has  illustrated  the  basic  simplicity  of  the  plan  and 
the  feasibility  of  programming  such  a  plan  for  digital  computation.; 
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Some  salient  features  are: 

a.  The  limit  on  size  and  complexity  of  the  system  to 
be  analyzed  is  that  caused  by  capacity  of  the  com- 
puter,  more  specifically,   the  capacity  of  the  partic- 
ular method  chosen. 

b.  The  time  required  for  analysis  is: 

1.  Read  load  deck  -  -  5  minutes 

2.  Computation  of  all  the  terms  of  the  performance 

equation.     This  will  vary  with  the  complexity  of 
the  system.     For  the  system  used  in  example 
number  1,    the  time  was  approximately  8  minutes. 

3.  Computation  and  punching  the  performance  function. 
Again  this  time  will  vary  with  the  system.     For 
example     number  1  the  time  was  approximately  4 
minutes.' 

4.  Computation  of  frequency  response  for  the  system 
of  example  1  -  -  slightly  under  2  minutes  per  fre- 
quency point. 

As  a  future  project  it  is  recommended  that: 

a.  The  program  be  refined  to  reduce  computation  time. 

b.  The  refined  program  be  altered  to  make  use  of  aux- 
iliary equipment  such  as  tape  to  lessen  present  phys- 
ical limitations. 

c.  With  greater  storage  capacity,    the  performance  func- 

tion routine  be  extended  to  enable  computation  of  the 
performance  function  as  a  ratio  of  polynomials  of  S. 

d.  Having  accomplished  (c),   an  investigation  be  made 
into  the  feasibility  of  extending  the  program  to  include 
factoring  the  polynomials  and  obtaining  a  time  res- 
ponse. 
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APPENDIX  A 


Program  Summary  No.    1,   L     Subroutine  (M0075  -  M0585) 

The  purpose  of  this  summary  is  to  explain  the  details  of  the 
program  covering  the  method  explained  in  paragraph  3.4.     Using 
the  G-array  of  paragraph  3.  3.  ,   the  first  summer  listed  is  located 
by  knowing  the  number  of  components  and  the  "width"  of  the  array. 
For  the  sample  system  of  paragraph  3.3.,   this  starting  point  is 
Gpn,   the  value  of  which  is  (6).     As  in  paragraph  3.4.  ,    the  G-array 
is  searched  by  rows,    starting  with  G,    (encountering  a  zero  indicates 
the  end  of  the  row)  until  the  search  number  is  found.     It  is  located  in 
Gr.     To  find  the  associated  output  (first  number  of  the  row  containing 
G5)  the  equation:     L  =  M  INTEGER  (L/ M  is  used;   M  is  the  width  of 
the  array  and  L.  is  the  G  subscript  (5).     A  new  value  of  L  is  computed 
then  as  (4)  hence,    the  number  wanted  is  G4  =  (2).     Continuing  as  in 
paragraph  3.4.,    no  more  outputs  of  (6)  are  found.     The  (2)  is  used  as 
a  search  number  and  is  found  to  go  to  (5)  and  also  negatively  back  to 
(7). 

Before  proceeding  further,   it  is  best  to  look  at  the  XX-array. 
(The  numbers  in  parenthesis  indicate  the  subscript  of  XX  locating 
the  position  in  the  array). 


XX-ARRAY 

(0) 

6 

(1)      (2) 

6         2 

(3)      (4) 

(5) 

6         2 

5 

(6)      (7) 

(8) 

-6         2 

7 
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Note;    that  the  (-)  appears  on  the  first  number  of  the  line  rather  than 
on  the  number  that  caused  it.     The  command  that  finds  the  search 
number  in  the  G-array  (M0230  and  M0240)  checks  for  both  positive 
and  negative  numbers  and  if  negative  sets  a  variable  R  =  1  (0255). 
Then  when  the  line  is  recorded,    each  term  is  multiplied  by  R  (M0350). 
After  the  first  number  of  the  line  is  computed,    R  is  set  to  +1.     This  is 
desirable  sihce  in  later  computations,    it  is  the  sign  of  a  loop  that  is 
important  ,     so  the  sign  of  the  firvst  number  only  is  checked. 

The  method  of  this  routine  is  adequately  described  in  paragraph 
3.4.     The  following  is  concerned  primarily  with  the  mechanics  of  the 
program.     As  noted  earlier,   before  recording  a  line  in  the  XX-array 
certain  checks  must  be  made,   the  first  of  which  is:     "Check  of  a  com- 
plete loop  (M02  65)".     This  checks  the  next  number  to  be  recorded 
against  the  first  number  of  the  line  and  if  they  are  the  same,    the  num- 
bers of  that  line  represent  units  comprising  a  complete  loop  and  are 
placed  in  the  ZZ  list  (M0500).     The  next  check  is  to  see  if  the  point 
is  already  recorded  in  the  appropriate  line.     The  newly  found  number 
must  be  checked  against  all  numbers  of  the  line  (M0275  to  M0285). 
If  it  is  a  duplicate  point,    a  new  line  is  not  recorded  but  the  search  is 
continued.     If  not  a  duplicate  point,    then  it  must  be  checked  against 
all  summers  previously  used  as  starting  points;   these  are  recorded 
by  (M0485)  and  the  check  is  made  by  (M0320). 

The  mechanics  of  recording  lines  in  the  XX-array  is  fairly  com- 
plicated and  is  most  easily  understood  after  getting  farther  along  in  the 
array,   although  the  same  commands  are  used  to  record  the  first  few 
lines.     Consider  the  stage  of  operation  in  paragraph  3.4.  ,    at  which  is 
recorded: 


line  (e)  -6273  (XXg  to  XX      ) 


(9) 
-6 

(10) 
2 

(11) 

7 

(12) 
3 

(13) 

(14) 

(15) 

(16) 

-6 

2 

7 

4 

line  (f)  -6  2  7  4  (xx13  to  XX16) 

(17)      (18)      (19)      (20) 
line  (g)  -6  2  5  8  (XX1?  to  XX2Q) 
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(2  1)      (22)      (23)      (24)      (25)  (XX     to  XX      ) 

line  (h)  -62581 

(4)  (26)      (27)      (28)      (29)      (30) 

line  (i)  -6  2  7  4  8  (XX2gtoXX3()) 

(31)      (32)      (33)      (34)      (35)    (36) 

(5)  line  (j)  -6  2  7  4  8  1       (XX^toXXgg) 

line  (k) 

The  array  is  separated  into  "blocks"  by  the  solid  lines  and  the  blocks 
are  given  numbers  indicated  by  the  number  to  the  left  of  the  array. 

To  handle  this  recording  procedure,    the  following  variables 
are  used: 

W  Number  of  current  block  (5) 

Y  The  XX  subscript  of  the  first  number  of  the 
last  line  of  the  previous  block 

V  The  XX  subscript  of  the  first  number  of  a 
line  containing  the  present  search  number. 

FF  Present  search  number 

N  In  general,    the  XX  subscript  of  the  first 

number  of  the  last  recorded  line. 

From  the  array  above,    assume  line  (j)  has  just  been  recorded  using 
FF  =  XX„n  =  (8)  as  a  search  number;   then 


w 

=  5 

V 

=  26 

N 

=  31 

Y4 

=  17 

Yr 

=  26 

Returning  to  search,    no  more  outputs  of  (8)  are  found  and  it  is  desired 
to  pick  up  XX25  =  (1)  as  the  new  value  of  FF. 


M0400  V  =  V  -  W  =  21 


FF  =  XXy  +  w  _  1  =  XX25  -  (1) 
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Return  to  search  the  G-array  and  find  (-1)  going  to  (7)  and  after 
making  the  3  checks  noted  previously  increase  NtoN  +  W  +  l  =  37 
and  record  (M0350)  line  (k)  as 


+6,     2,     5, 


1,     7    (XX3?  to  XX42) 


No  more  outputs  from  (1)  are  found  and  a  new  search  number  must 
be  picked  up.     From  the  array  above,    it  is  seen  that  it  is  not  desired 
to  go  back  another  line  but  to  start  a  new  block.     The  reason  for  the 
Y,    v  variable  is  now  apparent.     Change  V,   as  before,    toV=V-W=16 
and  check  it  against  Yw    1   to  see  if  block  4  has  been  finished  (M0405). 


V  =  N  =  37 
FF  =  XX 

W  =  W  H 


-  "X"X~         -   7 
V  +  W  42 


Y6=37 

N  =  N  -  1 


(The  reason  for  this  is  not  apparent  at 
present  but  will  be  explained  shortly) 


The  variables  are  now  set  up  to  go  back  to  search  and  complete  anr 
other  block.     The  process  continues  until  the  XX-array  appears  as 


(4) 


(5) 


(6) 


line  (h) 
line  (i) 


line  (j) 
line  (k) 


line  (1) 
line  (m) 


-62581 
-62748 


-627481 
+625817 


6    2    5     8    17    3 
6    2    5    8    17    4 


(XX 
(XX 

(XX 
(XX 

(XX 
(XX 


21-25 
26r-30 


31-36 
37-42 


43-49 
50-56 


Note  that  nothing  from 
line  (j)  is  recorded  here 
since  (1)  goes  to  (6)  and 
forms  a  complete  loop; 
also  goes  to  (7)  which  is 
a  duplicate  point. 
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Having  recorded  line  (m)  and  set  up  to  start  a  new  block,    the  values 
of  the  variables  are: 

V  =  N  =  50 
FF  =  4 

W  =  W  +  1  =  7 
Y?  =  50;   Y6  =  37 
N  =  N  -  1  =  49 

Return  to  search  with  FF  =  4  and  find  it  goes  only  to  (8),   a  duplicate 
point. 

V  =  V  -  W  =  43 

FF  =  XXy  +  w  _  j  =  3 

Return  to  search  and  find  (3)  goes  only  to  (6),   a  complete  loop,    so 
nothing  is  recorded. 

V  =  V  -  W  =  36 

V  is  less  than  Yg  hence  it  is  desired  to  start  a  new  block 
but  the  routine,   for  this  summer,   is  finished  since  there  is  nothing 
recorded  in  block  (7).     Recall  that  if  anything  had  been  recorded  in 
block  (7),    N  would  have  be^n  changed  to  N  +  W  +  1  =57.     Now, when 

V  is  set  equal  to  N  =  49,    the  check  against  Yw  indicates  the  end. 
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APPENDIX  C 


The  numbers  on  the  preceding  pages  of  this  Appendix,    represent 
the  component  numbers  of  each  term  of  the  performance  function  of 
the  sample  system  of  Chapter  4.     To  illustrate  the  complexity  of  the 
system  more  clearly,    the  performance  function  is  written  in  the 
more  conventional  form  below.     Since  there  are  so  many  terms  the 
numerator  and  denominator  are  written  separately. 

NUMERATOR: 

G6  "  G6  G2  +  G6  G3  +  G6  G4  +  G6  G4  G5  "  G6  G2  G3  "  G6  G2  G4  " 

-G6  G2  G4  G5  +  G6  G3  G4  +  Gg  G3  G4  G5  -  Gg  G2  G3  G4  -  Gg  G2  G3  G4  G5 
+  G?  G3  G4  G5  Gg  -  G?  G3  G4  G5  Gg  G2  +  Gt  G2  G3  G4  G5  Gg 

DENOMINATOR : 

1  +  Gg  +  G7  G3  G4  G5  Gg  +  Gx   G2  G3  G4  G5  Gg  -  G2  +  G3  +  G4  +  G4  G&  + 

+  Gg  -  Gg  G2  +  Gg  G3  +  Gg  G4  +  Gg  G4  G5  -  G?  G3  G4  G5  Gg  G2  - 

"  G2  G3  "  G2  G4  "  G2  G4  G5  "  G2  G6  +  G3  G4  +  G3  G4  G5  +  G3  G6  + 

+  G4  Gg  +  G4  G5  Gg  -  Gg  G2  G3  -  Gg  G2  G4  -  Gg  G2  G4  G5  +  Gg  G3  G4  + 

+  Gg  G3  G4  G5  -  G2  G3  G4  -  G2  G3  G4  G&  -  G2  Gg  Gg  -  G2  G4  Gg  - 

-  G2  G4  G5  Gg  +  G3  G4  Gg  +  G3  G4  G5  Gg  -  Gg  G2  G3  G4  -  Gg  G2  G3  G4  G& 

"  G2  G3  G4  G6  '  G2  G3  G4  G5  G6 
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